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^N) ' Abstract. In this paper, we are concerned with the system of the compress- 

ible Navier-Stokes equations coupled with the Maxwell equations through the 
Lorentz force in three space dimensions. The asymptotic stability of the steady 
state with the strictly positive constant density and the vanishing velocity and 
electromagnetic field is established under small initial perturbations in regular 
Sobolev space. For that, the dissipative structure of this hyperbolic-parabolic 
system is studied to include the effect of the electromagnetic field into the 
viscous fluid and turns out to be more complicated than that in the simpler 

P^ ' compressible Navier-Stokes system. Moreover, the detailed analysis of the 

^A ' Green's function to the linearized system is made with applications to derive 

^^ the rate of the solution converging to the steady state. 
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1. Introduction 



Plasma dynamics is a field of studying flow problems of electrically conducting 
fluids, particularly of ionized gases j321 I52j . The scope of plasma dynamics is very 
broad. A complete analysis in this field consists of the study of the gasdynamic field, 
the electromagnetic field and the radiation field simultaneously |50l [53] . Little has 
been done on such a complete analysis. Thus, many studies at both the microscopic 
kinetic level and the macroscopic continuous level are made in the individual field 
or in a little complicated situation taking into account the interaction of two fields 
with the rest one as a subsidiary factor [H |6l |32j [50l [52j [53] . 

At the macroscopic level, the motion of plasma fluids, for instance ionized gases, 
is described in terms of the classical gasdynamic quantities, such as density, pres- 
sure, flow velocity, etc, under the influence of the electromagnetic fleld induced by 
the electrically conducting fluid itself. The fundamental system to govern the time 
evolution of those macroscopic quantities consists of the gasdynamic equations, e.g. 
the Navier-Stokes equations in the viscous case, coupled with the equations of the 
self-consistent electromagnetic field, e.g. the Maxwell equations. In general, plasma 
is a mixture of various species: charged ions and electrons, and neutral particles. In 
many simplified cases when the variation of one composition is small and its effect 
is negligible, the plasma may be considered as a single fiuid. 

In this paper, we are interested in the motion of one fluid described by the 
compressible Navier-Stokes equations coupled with the Maxwell equations through 
the Lorentz force. In the iscntropic case, the motion equations read as 



(1.1) 



' 9t?i + V • (nu) = 0, 
dt{nu) + W ■ (nu (g) u) + VP{n) = -n{E + u x B) + I'Au, 
dtE ~V X B — nu, 
dtB + V X E = 0, 
W ■ E = lib - n, W ■ B = 0. 



Here, the unknowns are n = n{t,x) > 0, u = u{t,x) £ M.^, E = E{t,x) E M.^ and 
B = B{t,x) G M^ over {t > 0,x G R^}, denoting the fluid density, fluid velocity, 
electric field and magnetic field, respectively, j/ > is a constant denoting viscosity 
coefficient, and rib > is a constant denoting the uniform background density 
(e.g. of ions). P depending only on n denotes the pressure function with the usual 
assumption that P is smooth in the argument and P' > 0. Initial data of the 
system is given as 

(1.2) n(0, x) = no(x), w(0, x) = uo{x),E{0, x) = Eo{x), B(0, x) = Bo{x),x e R^ 

For convenience of presentation, we call (|l.ip the Navier-Stokes-Maxwell system. 
Notice that the same terminology was used in Masmoudi [47] but for a different 
modelling system. 

The consideration of system (jl.ip is motivated by three issues that we shall 
address in more detail later on. The first issue is to generalize some known results 
of the viscous fiuid, cf. Matsumura-Nishida [5l], Kawashima [55], Hoff-Zumbrun 
[T7] . to include the effect of the electromagnetic field. Notice that system (|1.1|) 
is of the hyperbolic-parabolic type. Indeed, from the later proof, the nonlinear 
coupling of system (jl.ip makes it possible for the hyperbolic Maxwell system to 
gain some weak spatial regularity or energy dissipation rate from the diffusion 
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effect of tlie viscous fiuid. On tfie other liand, as far as the large time behavior is 
concerned, the speed of convergence for the fluid velocity has to be reduced because 
of the appearance of the electromagnetic field. The second issue is to compare 
the current study of the Navier-Stokes-Maxwell system with that for other two 
relative models, the Navier-Stokes-Poisson system and the magnetohydrodynamic 
equations, cf. Li-Matsumura-Zhang [35], Umeda-Kawashima-Shizuta I61j . The 
final issue is to expose the special dissipative structure and regularity-loss property 
of system (|l.ip . which have been recently investigated for both fluid and kinetic 
models, cf. Hosono-Kawashima [121, Ide-Haramoto-Kawashima [52], Duan [TU] . 
Duan-Strain [13], Duan [9]. Notice that we shall derive system (|1.1|) from the 
one-species Vlasov-Maxwell-Boltzmann system [^ by employing the Liu- Yang- Yu's 
macro- micro decomposition |38| . This also provides a link between the macroscopic 
Navier-Stokcs-Maxwell system and the kinetic Vlasov-Maxwell-Boltzmann system 
for the single fluid with the electromagnetic held; the similar dissipative property 
for both systems will be also discussed later on. 

1.1. Main results. The Navier-Stokes-Maxwell system (|l.ip has a class of natu- 
rally existing constant steady states in which the density is Tib, the magnetic fleld 
is an arbitrary constant vector B^ and both the velocity and the electric field are 
zero. Here and hereafter, we call Bb the uniform background magnetic fleld. In 
this paper, B\^ = is assumed, and we are concerned with the asymptotic sta- 
bility of the constant steady state with the density rib and the vanishing velocity 
and electromagnetic field for the Cauchy problem (|l.ip . (|1.2p . We shall point out 
later on that part of results can be directly generalized to the case of non- vanishing 
constant magnetic fleld, that is, B\, ^ 0. 
To the end, set 



7 = V-P'(nb) > 0, ^ = VH^ > 0, /i = i//?ib > 0, 
and take change of variables 

Then, the Cauchy problem (|l.ip . (|1.2p can be reformulated as 



(1.3) < 


dt 


/ 
V + 7Vp + I3E - 


^V • {pv), 

- /iAu = ~-^v ■ Vw ( Vp) 

7 ~ 11 

— -V X B + v{ )Av, 

P P + Uh Uh 


and 








dtE -V X B - I3v ^ -pv, 

dtB + S/xE = 0, 

V-E = --p, V-i? = 0, 
I 7 


(1.4) 








with 










(1.5) ^ 


(0, 


x) = po{x), 


z;(0. 


x) = vo{x),E{0,x) = Eo{x),B{0,x) = Boix),x € R\ 
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One of the main results concerning the global existence and large time behavior 
of solutions to the above Cauchy problem on the Navier-Stokes-Maxwell system is 
stated as follows. 

Theorem 1.1. Let N > A. Assume that initial data [pQ, vq, Eq, Bq] satisfies the last 
equation of (jl.4p and || [po, ^o, ^Oj ^o]||ff" ^s sufficiently small. Then the Cauchy 
problem (|1.3p . (|1.4I) . (|1.5p admits a unique global solution [p,v,E^B] with 

[p,z;J,B]GC([0,oo);i7^(M3)), 

p e L^iO, oo); iJ^(R3)), Vv £ 1^(0, oo); H^iR^)), 



(1.6) \\[p{t),v{t),E{t),B(t)]\\j,.+ {\\pis)\\j,. + \\Vv{s)\\j,. 

Jo 

+ ||V^(.s)||^„_, + \\\/^B{s)\\l^^s)ds < C\\[po,vo,Eo,BoWh. 

for any t > 0. Moreover, if \\[po,vo, Eq, Bq]\\i1(-^jjn+2 is further sufficiently small, 
then the obtained solution [p,v, E, B] satisfies 



Mm 


<(l+i) 


Mm 


<(l+i) 


mm 


<(l + i) 


mm 


<(l+i) 



3ln(3 + t), 



for any t > 0. 



We now give some remarks on the above theorem. First of all, the global ex- 
istence result stated in Theorem 11.11 also holds for the case oi B\, ^ 0. In fact, 
the proof is easily modified to derive the a priori estimate as in (|1.6p by noticing 
that the additional linear term generated by i3b 7^ is m x i?b which is orthogonal 
to the velocity u; see the proof of Theorem 16.11 Secondly, (jl.6p shows that the 
electromagnetic field satisfying the Maxwell system indeed has some time-space 
integrability property weaker than the fluid components, whereas its highest-order 
spatial derivative is not time-space integrable, which becomes a typical feature of 
this kind of systems with the regularity-loss property. 

Thirdly, the time decay rates stated in Theorem 11.11 depend essentially on the 
analysis of the Green's function G of the linearized Navier-Stokes-Maxwell system; 
see Theorem 13. 1[ Theorem 15.11 Theorem 15.21 and Theorem 16.21 Specifically, from 
Theorem 15.11 the Fourier transform of G have different behaviors over different 
frequency domains, in the rough way as 

\k\<e, 

(1.7) \Git,k)\<{e~'^^'>* e<\k\<L, 

(i)t 

^ |A:|>L, 

for two properly chosen constants 0<e^l<Ci<CK); more precise estimates on 
G can be found in Theorem l5.11 It should be emphasized those time rates obtained 





'e-0(i)|tf t +£-0(1)1/01"* 


\Git,k)\<< 


e-0(i)t 




g-0(l)t ^g-0(l)|fc|^t_^ 
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in Theorem 12.21 on the basis of the elementary Lyapunov property of the linearized 
system, which implies 



til 

JTIT 



(1.8) \Git,k)\<e (i+i'=i^)% fceM^ 

fail to yield the ones of the nonlinear solution given in Theorem 11.11 above . Notice 
that (jl.Sp . which is a simpler form of (jl.7p . only provides the partial information 
of G. 

Finally we point out some differences of time rates between the linearized and 
nonlinear cases. From Corollarv 15.11 the solution [p,v,E,B] to the linearized ho- 
mogeneous system corresponding to (|1.3p - (|1.4p decays as 



Mm 


<(i+t) 


Mm 


<(l+i) 


imm 


<(l+i) 


\\B{m\ 


<(l + i) 



for any t > 0, provided that [po,vo,Eo,Bo] belongs to L^ (1 H^ for properly large 
s. Thus, in Theorem 1 1.1[ the time decay rates of v and B are optimal in the sense 
that they are the same as those in the linearized case, but other two rates for p 
and E are not optimal because they are strictly smaller than those in the linearized 
case. From the proof of Theorem 1 1 . 1 1 later on, the slower time decay for p and E 
mainly results from the nonlinear effect, particularly in the presence of a nonlinear 
term v x B in the momentum equation of system (jl.3p . 



1.2. Motivations and literature. As mentioned before, there are three motiva- 
tions resulting the current investigation of the Navier-Stokes-Maxwell system. The 
first motivation is to generalize some known results of the viscous fluid, particularly 
about the global existence and large time behavior of classical solutions near non- 
vacuum steady state, in order to include the effects of the electromagnetic field. 
In fact, there are extensive studies concerning those results for the Navier-Stokes 
equations in the context of gasdynamics, in the form of 

r dtn + V • {nu) = 0, 

1 dt{nu) + V • {7iu (g) u) + VP{ti) = i^Au, 

where the electromagnetic field is absent and thus the Maxwell system is decoupled 
from the fluid dynamic equations. Wc here mention only some of results related to 
our interest. 

One of the most important aspects of the Navier-Stokes system (|1.9p is due to its 
hyperbolic-parabolic character. It is crucial to understand the interaction between 
the hyperbolic conservative part and the parabolic diffusive part. Indeed, as we 
shall see even in the general situation ^E\ [S71 [SS], this kind of interaction can 
often induce some weak smoothness and dissipation of the hyperbolic component 
of the coupled system. It is Kanel who first observed this property in the proof 
of global existence of solutions to p.9p in one space dimension. Later on, in three 
space dimensions, global well-posedness of the near-constant-equilibrium solutions 
to the Cauchy problem or the initial boundary value problem was established by 
Matsumura-Nishida [48l [49l [51] using the classical energy method. Since then, the 
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perturbation theory of the Navier-Stokes system has been investigated by many 
authors. 

The study of the dissipative property of the system (jl.91) is not only important 
for the global existence of solutions under small perturbations but also useful for 
obtaining the rates of convergence of solutions trending towards the equilibrium. 
Still in [49], Matsumura-Nishida obtained the rate of convergence for the Navier- 
Stokes system in K^ as 

\\[n-n^,u]\\H2 <C(l+t)-^ 

if initial perturbation is small belonging to H'^OL^. Ponce [S3] gave the time decay 
rate in L^ with 2 < p < oo as 

\\[n-m„u]\\LP<C{l + t)-"^^'-l^^ 

provided that small initial perturbation belongs to H^ D W^'^ for s large enough. 

On the other hand, a general approach for obtaining the optimal time decay 
of solutions in LP space with p > 2 in any space dimension was developed by 
Kawashima [28] [26] and Shizuta-Kawashima [57] . Precisely, the main tool used by 
Kawashima [281 126] is the Fourier analysis applied to the linearized homogeneous 
system, and the key part in the proof is to construct some compensation function 
to capture the dissipation of the hyperbolic component in the solution and then 
obtain an estimate on the Fourier transform U of the solution U as 

(1.10) \U{t,k)\<Ce~^^'\Uo\ 

for any t and k. Thus, from the above estimate, it is well known the solution over the 
high frequency domain decays exponentially while over the low frequency domain it 
decays polynomially with the rate of the heat kernel. Moreover, in [57], the authors 
proposed a condition, which is now called the Shizuta-Kawashima condition, in 
order to assure the time decay of solutions to linear systems of equations of the 
hyperbolic-parabolic type. The Shizuta-Kawashima condition also plays a key role 
in the study of stability and large time behavior of the solution to the nonlinear 
system; see Hanouzet-Natalini [21], Yong [66], Bianchini-Hanouzet-Natalini [3] and 
references therein. For some extensions of the Shizuta-Kawashima condition, refer 
to Beauchard-Zuazua [2] and Mascia-Natalini [46] . 

About the time decay in L^ space with 1 < p < 2, Zeng [67] provided a com- 
plete analysis of the Green matrix with some sharp, pointwise bounds for the one- 
dimensional viscous heat-conductive fluid system. The difficulty for the rate in LP 
for p close to one lies in the fact that the density satisfies the mass conservation 
law. The result of [ST] was later extended by Liu-Zeng [32] to the case of zero heat 
conductivity and also by Liu-Zeng [33] |33] to general hyperbolic-parabolic systems 
of conservation laws in one space dimension. Besides, an extension of j67) to the 
case of multi-dimensions was made by Hoff-Zumbrun [17l[18]. Precisely, [17] showed 
that for the Navier-Stokes system in R'^, 

( Cri^'^~p\ 2<p<oo, 
\\[n- nb,u] II LP < < 

[ct-^'^^-p\ l<p<2, 

for all large t > 0, provided that the small initial perturbation belongs to H'' n L^ 
with s > 4. The main idea of [T7] is to write the Green matrix G in terms of the 
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Fourier transform as 

Git,k) = G^^\t,k) + G^^\t,k), 
&^\t, k) ^ G{t, k)x{k), G^^\t, k) = G{t, k)[l - x(fc)] 

with x(fc) a- smooth cutoff function supported on the finite domain, and tlius sliift 
the difhcuhy of the analysis to the high-frequency part G*-^-* (i, k). We also mention 
the work in the case of multi-dimensions by Liu- Wang |37| for the pointwise estimate 
on the Green's function, and similar work by Li |34] and Wang- Yang [64]. Finally, 
we mention a recent work Li-Zhang |36j for the study of the optimal large time 
behavior of the Navier-Stokes system in the Besov space when initial perturbation 
does not belong to i^ but some homogeneous Besov space which can induce some 
extra time rate. 

In addition, when there is a given potential force, the density of the naturally 
existing steady state is not a constant but a function of spatial variable. In this 
case, the Fourier analysis fails due to the difhculty of variant coefficients. The slow 
time decay property of solutions was studied by Deckelnick [3 [8] on the basis of the 
refined energy estimates together with an ordinary differential inequality. Notice 
that even though the time rate obtained in [71 [8] is not optimal, the developed 
method has proven very useful with applications to the case when the spectral 
analysis which can lead to the optimal time decay rate is not clear; see Yang-Zhao 
[55] for a study of the Vlasov-Poisson-Boltzmann system. An extension of [7] H] 
was made by Ukai- Yang- Zhao [60] to obtain the almost optimal time decay rates. 
Recently, Duan-Ukai- Yang-Zhao [15] developed a method of the combination of 
the spectral analysis and energy estimates to deal with the optimal time decay of 
solutions, and an extension of [15] was made in [11]; see also [M] for the detailed 
presentation of the method. When the given external force depending only on 
spatial variable has non-zero rotation, the velocity of the steady flow thus can not 
vanish; for such case, the existence and time asymptotic stability of the nontrivial 
steady state were studied by Shibata-Tanaka |55| and [56], respectively. 

At this moment, those mentioned works about the optimal time decay of solu- 
tions over the unbounded domain are concerned with the case of the whole space. 
Kagei-Kobayashi (25l [24] proved the large time behavior of solutions over the half 
space in M^. Kobayashi-Shibata [31] and later Kobayashi [30] gave an analysis 
of the Green's function and large time behavior of solutions for the compressible 
Navier-Stokes equations in an exterior domain of MP. 

The second motivation is to compare the current study of the Navier-Stokes- 
Maxwell system to the one of some modelling systems such as the Navier-Stokes- 
Poisson system in semiconductor theory and the MHD system in magnctohydro- 
dynamics. Here, we also mention some previous work of these two systems only 
related to our study. First of all, the Navier-Stokes-Poisson system can be used 
to simulate, for instance in semiconductor devices [45], the transport of charged 
particles under the influence of the electrostatic potential force field, and in the 
compressible case it takes the form of 

dtn + \/ ■ {nu) = 0, 

(1.11) { dt{nu) +V ■{nu®u)+ VP{n) = -nE + z^Am, 

S = V^, V-E = nb-n. 
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Here, compared with the Navier-Stokes-MaxweU system, the magnetic field B is 
omitted and the electric field E is generated by a consistent potential function $ 
therefore satisfying the Poisson equation 

A<1> = nb — n. 

It is easy to see that the compressible Navier-Stokes-Poisson system (|1.11[) is 
a coupled system of the hyperbolic-parabolic-elliptic type. This essential feature 
makes the Navicr-Stokcs-Poisson system have some special dissipative and time 
decay properties different from both the Navier-Stokes system and the Navier- 
Stokes-MaxweU system. Li-Matsumura-Zhang [3S] made a delicate analysis of the 
linearized Navier-Stokes-Poisson system and investigated the spectrum of the linear 
semigroup in terms of the decomposition of wave modes at lower frequency and 
higher frequency respectively, which can be used to deal with the influences of 
the rotating effect of the electric field. Precisely, for the rate of convergence, they 
showed that the solution over M^ satisfies 

||?i-nb|| <C{l + tri, 

\\u\\<C{l + tyi, 

\\E\\<C{l + t)-K 

for any t > 0, provided that initial perturbation belonging to H'^ (IL^ for s properly 
large is sufficiently small, and furthermore, these time rates are optimal due to the 
fact that they can also be obtained as the lower bounds if the Fourier transform of 
initial density perturbation riQ — rib is strictly positive near some neighborhood of 
the origin, i.e. 

(1.12) inf |n^r^b(0,fc)| >0. 

IfeKi 

As pointed out in [35j, it is the rotating effect of electric field which reduces the 
time decay rate of the momentum component compared with the Navier-Stokes 
system (jl.9p . We remark that if uq — nb(0,fc) behaves instead like |fc| as |fc| tends 
to zero, for instance 

(1 + |a;|)|no(a;) — Tn^\dx < oo 
i 

holds true, then those time decay rates obtained in j35j could be improved. More- 
over, ()1.12p implies that Eq ~ V^o does not belong to L^(E^), otherwise it is 
possible to improve the time rate of ||£'|| as (1 + i)^3. 

However, in our discussion of time decay rates for the linearized Navier-Stokes- 
MaxweU system, since the electric field E satisfies a time evolution equation and 
initial data Eq is given belonging to L^, the higher time rate (1 -|- f)~3 of \\E\\ is 
natural. In addition, the time rate (1 +t)~'s of ||m|| for the momentum component 
is still strictly slower than (1 -I- t)~* in the case of the Navier-Stokes system. This 
point essentially results from the magnetic field effect of B which has the slowest 
time decay rate among all the components of the system. Finally, ||n — n(,|| decays as 
(1 -|- t)~i in the fastest way, because on one hand, the constrict condition n — nb = 
—V • E makes n — ni, gain one more spatial derivative, and on the other hand, n 
at the linearized level turns out to be decoupled with and hence not effected by 
the magnetic field B with the slower decay, where the decoupling property is due 
to the hidden assumption in our study that the background magnetic field is zero. 
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Thus, it would be quite interesting to extend the current resuhs to the case of 
non- vanishing constant background magnetic field. 

After [35], Wang-Wu [B^ obtained the pointwise estimates of solutions to the 
Navier-Stokes-Poisson system through the analysis of the Green's function. Zhang- 
Li-Zhu [68] generalized the iscntropic case to the non-isentropic case with additional 
efforts taking care of the temperature equation. Hao-Li [16] proved the global 
existence of solutions in the Besov space. We also mention the work Zhang- Tan 
[69j for a study in two space dimensions and a survey Hsiao-Li |20) . 

Next, let us discuss a little about the magnetohydrodynamic system [11 [6]. In 
the isentropic case, it reads 

dtn + V ■ (nu) = 0, 

(1.13) { dt{nu) + V ■{nu®u) + VP{n) = (y x B) x B + vAu, 

dtB + V X {V X B - u X B) = 0, V ■ B = 0. 

Different from the Navier-Stokes-Poisson system, system (|1.13p models the in- 
teraction between the viscous fluid and the variable magnetic field. In one and 
two dimensions, Kawashima-Okada [29] and Kawashima [27], respectively, proved 
the global existence of smooth solutions near constant equilibrium states. Later, 
Umeda-Kawashima-Shizuta [HI] gave a detailed study of the linearized electro- 
magneto-fiuid system, where more general situation of non- vanishing constant back- 
ground magnetic field was considered. The authors obtained the same estimate on 
the solution as (|1.10p . which therefore implies the same time decay property as in 
the case of the Navier-Stokes system. Notice that the spectrum analysis was also 
made in [61]. We finally mention two recent works Chen- Tan [5] and Zhang-Zhao 
[70] both about the time decay rates for the nonlinear system, and Masmoudi |47] 
for the study of the global existence and exponential convergence rate of regular 
solutions to the full magnetohydrodynamic system over two dimensional bounded 
domain. Wc remark that the current method could be applied to the model in |47| 
to obtain some similar results in the framework of small perturbation. 

The third motivation is based on recent investigations on the fluid and kinetic 
systems in the presence of the electromagnetic field, such as the Euler-Maxwell sys- 
tem with relaxation [10] and the Vlasov-Maxwell-Boltzmann system [TS] [9] . For the 
linearized dynamics in all works }10[ I13[ [9] and the current work, the main idea is to 
construct the time-frequency Lyapunov functional to capture the dissipation of the 
electromagnetic field. Since the pure homogeneous Maxwell system preserves en- 
ergy and hence has no dissipation, the recovered dissipation of the electromagnetic 
field for those coupled systems is truly nontrivial. 

The optimal large time behavior of solutions near equilibriums for the Vlasov- 
Maxwell-Boltzmann system was proved by Duan-Strain [13] in the two-species case, 
but it is still unclear in the one-species case even though Duan j^ made the full 
linearized analysis of the one-species Vlasov-Maxwell-Boltzmann system and proved 
the global existence of classical solutions under small perturbation. In fact, the one- 
species Vlasov-Maxwell-Boltzmann system at the linearized level shows the same 
time- frequency Lyapunov property as in (|2.5I) . that is. 



dt£{U{t, k)) + c^^-^l^£{U{t, k)) < 0. 
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Thus, the hnearized solution operator also has the same time decay property as 
stated in Theorem 12.21 But, those rates are too slow to deduce the explicit time 
rate of solutions to the nonlinear system by using the usual bootstrap procedure. 

Different from the case of the one-species Vlasov-MaxwcU-Boltzmann system, 
our current study of the Navicr-Stokcs-Maxwcll system has made full use of the 
pointwise estimate on the Green's function of the linearized system. Our analysis 
of the Green's function here strictly refines the time decay property of all the 
components except the magnetic field, which is hence able to be used to obtain the 
time decay of the nonlinear solution. For the computation of the Green's function, 
we use the technique of [lOj to write the linearized system as two decoupled system, 
which makes it possible to separately solve them and also shift the whole difficulty to 
the analysis of some ordinary differential equations. On the other hand, the estimate 
on the Green's function in this paper is much more complicated than that in the 
case of the relaxed Euler-Maxwell system [TU] because of the hyperbolic-parabolic 
character of the Navier-Stokes-MaxwcU system, which results that multiple real 
eigenvalues could occur over the finite frequency domain. We shall make a careful 
analysis of the eigenvalues to overcome this difficulty. Therefore, the current method 
could also provide a clue to hopefully derive the optimal large time behavior of 
solutions to the one-species Vlasov-MaxwcU-Boltzmann system [3] , particularly the 
corresponding Green's function as studied by Liu-Yu J401I41| for the pure Boltzmann 
equation. 

Finally, it should be pointed out that all the coupled systems considered here 
and in [9l Uni E] have a common feature that they are of the regularity-loss type. 
A typical estimate on the solution U to the linearized system with this kind of 
regularity-loss property can be written as 

\Uit,k)\<Ce-*^''^'\Uo{k)\ 

for all t and fc, where (/>(fc) is a strictly positive, continuous and real- valued function 
satisfying 4>{k) — >■ as |fc| — >■ and oo both with the explicit polynomial rate in 
|fc|, see Lemma [2. II The regularity- loss of the system corresponds to the property 
of 0(fc) — >■ as |fc| — > oo, which is in general not true for the hyperbolic-parabolic 
system under the Shizuta-Kawashima condition. Therefore, the solution over the 
high frequency domain does not decay in time exponentially any longer but decays 
polynomially with any time rate depending on the regularity of initial data. Some 
recent studies of concrete models also have shown the similar phenomenon. Among 
them we mention the relaxed Euler-Maxwell system in plasma physics |59j . the 
Timoshenko system in biology [22l [23] and the hyperbolic-elliptic system in radia- 
tive hydrodynamics [19[ [551 \TI\ . Notice that [SI] considered the case of the nonzero 
constant background magnetic field. 

1.3. Notations and arrangement of the paper. Through this paper, C denotes 
a generic positive (generally large) constant, c denotes a generic positive (generally 
small) constant, and 0(1) denotes a generic strictly positive constant. For two 
quantities A and B, A < B denotes A < CB and A ^ B means cB < A < CB. 
For u = u{x), the Fourier transform u = u{k) of u is defined as 

u{k) = J^u{k) = / e~''''''^u{x)dx. 
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For an integer m > 0, ff™ denotes the Sobolev space H™{M.^) with the usual norm 
II • llffni. L^ = H'^ witli norm || • || is set when m = 0, and (•, •) denotes the inner 
product in L^. For a multi-index a ~ {ai,a2, as), we denote d" = d"^d"^d"^ , and 
the length of a is |a| = ai + a2 + cts- Finally, (V) is defined in terms of the Fourier 
transform as J"(V) = (C) := (1 + \£.\'^V^^. 

The rest of the paper is arranged as follows. In Section[21 we study the dynamics 
of the linearized homogeneous Navier-Stokes-Maxwell system by proving some Lya- 
punov inequalities, and use the Lyapunov property of the system in terms of the 
Fourier transform to obtain some elementary time decay estimates on the linearized 
solution operator. 

In Section 131 we explicitly solve the Green's function in the Fourier space. There, 
the key point is to analyze the equations of the electromagnetic field and shift the 
difficulty to solve a third-order ordinary equation. For the third-order ordinary 
equation, we study in Scction|4]some basic properties and the explicit representation 
of its roots, and obtain the asymptotic behavior of roots both at the origin and 
infinity. In Section [Sj we make a careful analysis of the pointwise bounds of the 
Green's function in the Fourier space, and use them to prove the refined time decay 
estimates on each component of the solution to the linearized system. 

In Section [6l we devote ourselves to the proof of the global existence of solutions 
to the nonlinear Cauchy problem, and apply the obtained linearized estimates to 
the nonlinear case in order to prove the large time behavior of solutions and hence 
complete the proof of the main result Thcorem ll.il 

In Appendix, we provide a derivation of the original Navier-Stokes-Maxwell sys- 
tem from the kinetic Vlasov-Maxwell-Boltzmann system by using the Liu- Yang- Yu's 
macro-micro decomposition. 

2. The linear homogeneous dynamics 
The linearized homogeneous system corresponding to (|1.3|) and (|1.4|) reads 
dtn + 7V ■ u = 0, 



(2.1) 



dtu + jVn + /3E - fiAu = 

dtB + V X E^O, 

V-E = --n, V-B = 0, i>0, xeM^ 

7 



Here and in the sequel, for simplicity, we use the same notations n, u, E, B as in 
(jl.ip to denote the unknown functions. In terms of the Fourier transform in x 
variable, the above system is equivalent with 

dtfi + "fik ■ u — 0, 





(2.2) 



dtu + jik 


h + [3E + fi\k 


^u 


dtE - ik 


xB' 


- Pu^O 




dtB + ik 
ik-E = - 


xE = 

— n, 

7 


= 0, 
k-B^ 


0, 



i > 0, fee 



1.3 



In this section, we make some elementary Fourier analysis as in |101 113[ [9] to exploit 
the dissipative structure of the system p.ip or cquivalently p.2|) by showing that 
the system has the Lyapunov property both in the Fourier and physical spaces. The 
main idea is to construct the time-frequency functionals and the induced instant 
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energy functionals that can capture the fuU dissipation of the system, particularly 
the dissipation of all the non-diffusion components in the solution. Moreover, basing 
on the obtained time-frequency Lyapunov inequality, we use the Hausdorff- Young 
inequality to study the time decay property of solutions in the standard way as 
in |28| . Those time decay estimates will be refined by further investigation of the 
Green's function of the linearized system in the next sections which are the main 
parts of this paper. 

2.1. Lyapunov property. First, this property is established with respect to the 
solution in the Fourier space in the following theorem. Throughout this subsection, 
we use (a | b) to denote the complex dot product of any two complex numbers or 
vectors a and b. 

Theorem 2.1. Let U := [n, u, E, B] satisfy the system ((2T|) Jor all t>0, x eR^. 
Then, there is a time-frequency functional £(■) such that 

(2.3) £{Uit, k)) - Pit, fc)|2 := \n{t, k)\^ + \u{t, k)\^ + \E{t, fc)|2 + \B{t, fc)|2 



(2.4) dtS{U{t, k)) + c\n{t, fc)p + c\k\^\u{t, k)\ 

\k? 



for allt>0,ke] 

(2.5) dt£{U{t,k)) + c 



(1 + |A:|2)2 
In particular, 



l^(t,fc)r-t-c 



ifci 



(l + |fc|2)3 



|B(t,fc)|2<0 



r£{Uit,k))<0 



'(l + |fc|2)3 

for all t> 0, k£ M.^ . 

Proof The Fourier transform representation (|2.2p of system (|2.ip can be written 



as 



/n\ 

u 
E 


+ 


1 -iik^ 

7ifc 




\ 


—iky. 




In 
u 
E 


\ 




\b) V ^fcx y \b) 


+ 


/O 00\ 
0/30 
0-/3 




/n\ 

u 

E 








\^ 





oy 




\b) 



/O 00\ fn\ 

0^i\k\^0 

00 
\0 00/ 



u 

E 

\B/ 



Since both the coefficient matrices of the first-order hyperbolic part and the zero- 
order relaxation part are skew-symmetric, the direct computation gives 



(2.6) 

From the identity 



-dt\[h,u,E,B]\''+^i\k\^\u\'' 



0. 



{ijkn I -fikn + I3E) = {-i^\k\'^ + I3'^)\n[' 
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and also using the second and first equations of p.2p to deduce 

{ijkh I jikh + f3E) = (ijkfi \ —dtu — ^\k\ u) 

= —dtiijkfi I u) + {ijkdtn \ ii) — 'yfj.{ikn \ \k\ u) 
= —dt{i^kh I -u) + 7 |fc • u| — 'yfi{ikh \ \k\ ii), 

one has 

dtiijkh I u) + (72|fc|2 + /32)|fi|2 = ^2|^ . ^|2 _ ^^(^ikh I Ifcpu), 

which after further taking the real part, using the Cauchy inequality for the last 
two terms and then dividing by 1 + |fcp, implies 

y{{ijkh I u) 
In the similar way, from the identity 



(2.7) a /T:::r +cN^<g|fcpH^ 



^2 

{E I jikn + /3E) = I3\E\^ + ^1^ ' -^1^ 



and also using the second and third equations of (|2.2p to derive 

(i? I -fikh + pE) = [E I -9tu - /i|fcpu) 

^-dt{E\u) + {dtE\u)~is\k\^[E\u) 

= -dt{E I It) + {ikxB + liu\u)- ^L\k\^{E I u), 

one has 

-72 . 
(2.8) at(-B I u) + /3|-Bp + -^|fc • ^|2 = /3|ft|2 + {ikxB\u)-- fi\k\''{E \ u). 

P 



Moreover, from the fourth and third equations of (|2.2p . 

|fc|2|B|2 = |fcxB|2 

= (ifc X B I ifc X S) 

= (ifc X B I a^i? - Pu) 

= dt{ik X B\E)- {ikx dtB \ E) ~ (3{ik x B \ u) 

= dt{ik X B\E) + \kx Ef ~ /3{ik x B\u), 

that is, 

(2.9) ~ dtiik X B \ E) + \k\^\B\^ = |fc x E\^ - l3{ik xB\u). 

We go on making estimates from the obtained identities (|2.8p and (|2.9p . By taking 
the real part of (|2.8p , using the Cauchy inequality for its right-hand last term and 
then multiplying it by |fcp/(l + |fcp)^, it follows 

(2-10) ^* (l + |fc|2)2 +^(i + |fc|2)2l^l 

^ /3|fcp 2 I g|fe|' 1^,2 , |A;p<H(zfc X B I u) 



(l+|fc|2)2' ' (1 + |A:|2)2' ' (l+|fc|2)2 

<C|fc|2|u 



2,„,2 \k\^m{ik X B \ ii) 



(l + |fc|2)2 
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On the other hand, from taking the real part of p.9p . using the Cauchy inequahty 
for its right-hand last term and then multiplying it by |fcp/(l + |fcp)'^, it follows 



(2.11) dt 



\k\'^m{~ikxB I E) 



< 



(l+|fc|2)3 
|fc|2 



(l + |fcP)3 



k X E[ 



•(l + |fcP)3 



\B\ 



C 



(l + |fcP)3 



\u\'< 



\k\' 



(i + lki^y 



\E\ 



c\km 



Define 



(2.12) £{U)^\[n,u,E,B]\^ 

d\{i"fkn I ii) 



■ Kl- 



1 



|fc|2fR(S|u) \k\^n{-ik X B \ E) 



(l + |fc|2)2 ' ' (1 + |A:|2)3 

with < Ki, K2 ^ 1 to be determined so as to guarantee that £{■) is the desired 
time- frequency functional satisfying p.3p and (|2.4p . Notice that as long as < 

Kl,K2 ^ 1, 

f(f/)^|[n,^,£;,B]|2, 



and hence (|2.3I) holds true. One can choose K2 > further small enough such that 
the sum of (|2.10|) and (|2.1ip multiplied by K2 gives 



(2.13) dt 



\k\^m{E I m) |fc|2fR(-ifc xB\E) 



(l + |fc|2) 



2\2 



K2- 



CK2 



(l + |fc|2)3 



+ (c- K2) 



iii;i- 



(l + |fc|2) 



2\3 



^2^ |/C|2in(jfc X i? I W 



(1 + l^p) 



(l + |fc|2)2' 

c\km 



2\2 



Since 

|A:|2«n(iA: x B \ u) 



< e 



\k\' 



(l + |fc|2) 



2^3 



ISI 



1 



|fc| = 



4el-h|/fc|2' 
< e- 



isr 



(l-f |fc|2)3 

holds for any e > by the Cauchy inequality, it follows from p.l3|) that 



i\k\'\u\' 



(2.14) dt 



|fc|2in(£; I ii) |fc|2<n(-ifc xi3\E) 

' K2 



(1-K|A:|2)2 

+ (c- K2) 



(l + |fc|2)3 
\k\^ 



ii + \k\^y 



\E\ 



--(irfrF'^'^^S'^'^'^'^- 



We now let K2 be fixed. By multiplying (|2.7p . (|2.14p by ki > further chosen 
small enough and then adding them to (|2.6p . (|2.4p follows. (|2.5p is the immediate 
consequence of (|2.4p . The proof of Theorem 12. II is complete. D 

By taking integration over k G K'^ for those time-frequency inequalities in Theo- 
rem [231 it is straightforward to obtain the corresponding estimates on the solution 
in the physical space, which describe the evolution of instant energy functionals 
with the optimal energy dissipation rates. Precisely, we have 
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Corollary 2.1. Let U := [n, u, E, B] satisfy the system (|2.ip for allt >0, x eR^. 
Let the time- frequency functional £{■) be obtained in Theorem \2.1\ Then, for any 
integer j > 0, 

^ / £{\k\^Uit,k))dk + c\\Vhif + c\\V^+'uf 

at Jjja 

+ c\\W^+^Wr^Ef+c\\W^+'^{Wr^Bf <0 
for any t > 0, where 

f £{\k\^U{t,k))dk^\\V^U{t)r. 
Particularly, for any integer N > 3, 



d_ ^ 
dt 



y f £{\k\^U{t,k))dk + c\\n\\j,^, + c\\Vu\\l, 



c1|v£;||^„-2 + c||v2b||^«_3 <o 



for any t>Q, where 

N 



J2 f £m'uit,k))dkr^\\u{t)rHr.. 



3=0 

Theorem 12.11 and Corollary 12.11 exactly show the dissipative structure of the 
Imearized Navier-Stokes-Maxwell system (|2.ip . The importance of this structure is 
that it exposes the dissipative property of the hyperbolic Maxwell equations in the 
full system, and also it plays a key role in the study of the nonlinear asymptotic 
stability of the constant steady state under small perturbations; see Theorem 16.11 
and its proof. 

2.2. Time-decay property. As in [TUl [T31 IS], it is now a standard procedure 
to derive from Theorem 12.11 the L^-L'' time decay property of solutions to the 
linearized system (|2.ip . Here, we consider it by obtaining the following lemma 
under a more general situation as in [3] , and also give the proof of the lemma both 
for completeness and for the later use in the proof of Theorcm l5.2l To do that, for 
f e M and 1 < s < 2 < (? < cx), define 

fO if ^ + 3(i-i) <0, 

^ s q ' ' 



m{£,s,q) 



if ^ + 3(i-i)>0,,s = g = 2. 



and i is integer, 

[£ + 3(i-i)] + l otherwise, 

where [•] means the integer part of the nonncgative argument. 

Lemma 2.1. Assume that for any initial data Uq, the linear homogeneous solution 
U{t) ~ A{t)Uo obeys the pointwise estimate 

\U{t,k)\<Ce-^^''^'\Uom 
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for all t > 0, k G M!^ , where (j){k) is a strictly positive, continuous and real-valued 
function over fc € R"^ and satisfies 



4>{k) 



0(l)|fc|'"+ as \k\ -^ 0, 
0(l)|fc|-'"- as\k\^oo, 



for two constants cr_ > cr+ > 0. Let j > be an integer, l<p,s<2<q<oo and 
£ > 0. Then, U{t) ~ A{t)Uo obeys the time-decay estimate 

(2.15) ||V^[/(t)|U, <C(l + t)"^'^"'^"^|j[/o|Up 

for any t > 0. 

Proof. Take a constant L > 0. From the assumptions on (t>{k), it is easy to see 

fc|A:r+ if|A:|<L, 



Take 2 < g < oo and an integer j > 0. By Hausdorff- Young inequality, 



(2.16) \\ViUit)\\L.<C 



J- -0(fc)t 



<C 



|fcpe-A|/cr+t^^ 



Ifc^e 



Li'{\k\<L) 



Un 



Li' 



Li'{\k\>L) 



h+h 



where - + 4- = 1- For /i, take 1 < p < 2 and by using the Holder inequality for 
^ = ^2^ + -L with p' given by i + ^ = 1, 

qi p' q/ p/ r- o •■' p p' ' 



h<c 



|A;Pe-^l*^l"+* 



It/n 



Pli^ ll-0||Lp'(|fc|<L)- 

LP -1 (|fe|<L) 



Here, the right-hand first factor can be estimated in a standard way |28j as 

_3 ri_i\_j_ 



\k\3^-MH'+t 



lV^ {\k\<L) 



< C(l + i) "+^'' '■'" '+ 



by using change of variable fci "+ — > fc, and the right-hand second factor is estimated 
by Hausdorff- Young inequalities as 

ll^o!lLp'(|fc|<L) < \\Uo\\lp' < C\\Uo\\lp- 

Therefore, for /i, one has 



Ii<C{l + t) '^^- '^ '^\\U4lp. 



To estimate I2 , take a constant ^ > so that 



l2=C 



|fcP'e^^l'=l"""*f7o 



<C sup |/fc|-^e-"l'^'l ""* 



|fcP'+'C/o 



L?'(|fc|>L) 



Here, the right-hand first factor decays in time as 

sup |/fc|-^e-'=l'=l"""* <C(l+t)" 

|fe|>L 
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We estimate the right-hand second factor as follows. Take 1 < s < 2 with - + ^ ^ I 

^ — — s s 

and take a constant e > small enough. Then, similarly as before, from Holder 
inequality for ^ = ^-—7- + -^, one has 

^ J q' g' q' g' 7 



IW+'Uo 



< 



L<!'(|fe|>L) 
|fc|-3(l+.) 



L''-i' (|fc|>i) 



,j+i+3{l+e}i 



^Un 



L^'i\k\>L) 



<a 



\k\ 



J + [«+3(i-i)] + 



Uo 



L-'(\k\>L) 



Here, due to s' > 2. further by Hausdorff- Young inequality. 



,u+ie+3(i 



')h 



Uo 



< 



L^'{\k\>L} 



|fcp + [^+3(i-i)]+^^ 



L=' 



Then, it follows that 



\k\'+'Uo 



Thus, I2 is estimated by 



Li'i\k\>L) 



< cj|v; 



<C^W^''^''--^''^Uoh^ 



i+[^+3(|-i)] + 



Uo\ 



Now, p.isp follows by plugging the estimates of /i and I2 into ()2.16p . This com- 
pletes the proof of Lemma 12.11 D 



By applying Theorem 12.11 together with Lemma 12.11 above to the linearized sys- 
tem (12.11). we have 



Theorem 2.2. Let 1 <p,s<2<q< 00, ^ > 0, and let j > be an integer. 
Let U = [n,u, E, B] satisfy the system (|2.ip for all t > 0, x gM.'^ with initial data 
U\t=o = Uo- Then, the solution of the linearized homogeneous system satisfies 

||V-''C/(i)||L, <C(l + i)-5(i-i)-i||[/o||^, + C(l + t)-i||V'"(^+^'^'«)t/o||L= 
for any i > 0. 



For the linearized system (|2.ip . Theorem 12.21 describes the unified time decay 
property of the full solution. But, it can not be directly applied to the nonlinear 
system (|1.3p - (jl.4p to obtain the time decay rates of the solution, particularly those 
rates for each component of the solution as stated in Theorem 11.11 Therefore, we 
turn to the study of the Green's function of the linearized system. 

3. Green's function 

In this section we explicitly solve the Cauchy problem on the linearized Navier- 
Stokes-Maxwell system p.ip or equivalently its Fourier transform (j2.2|l with initial 
data 

U\t^o^Uo = [no,uo,Eo,Bo], x eR^ 
which satisfies the last equation of (|2.ip . In fact, as in [10], p.ip can be written 
as two decoupled subsystems which govern the time evolution of n, V • m, SJ ■ E 
and V X u, V X i?, V X _B, respectively. For simplicity, we call two subsystems 
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by two systems for the fluid part and the electromagnetic part, respectively, even 
though only n and B can be fully decoupled. By solving some ordinary differential 
equations, we shall deduce the representation of the Green's function both for the 
fluid part and the electromagnetic part. 

3.1. Green's function of the fluid part. Wc flrst solve n. From (|2.2[) . by taking 
the time derivative of the first equation and using the second equation to replace 
dtii, one has 

dfn + -fik ■ {-jikn - jiE - /i|A:pu) = 0. 

Further by using ik ■ E ~ —-n and dtfi + jik • u = 0, the above is equivalent with 

dffi + fi\k\'^dtn + {-f^\k\^ + P^)n = 0. 
Thus, one has to solve the second-order ODE 

dffi + ^i\k\'^dtn + (72|fc|2 + ^2)^ ^ 0, 

(3.1) <( n(0,fc)=no(A:), 

^th[Q^ k) = —i^k ■ Uo{k). 
Consider the indicate equation 

(3.2) z^+n\k\''z + {-f^\k\^+/3^)=0, 
and its roots are denoted by 

1 ,.,o , 1 



(3.3) A±=A±(fc):=--M|A:p±-yMW^4(7N^+^. 



Set the representation of the solution h to (|3.ip as 

n(i, k) ^ C_ {k)e^- (*='* + C+ (fc)e^+ ('=)* 

except a finite number of values of |fc| where A_|_ = A_ become a two-multiple real 
root. Then, the initial conditions of p.ip imply 

rC_(fc) + C+(fc) = no(fc), 

1 A_(fc)C_(fc) + A+(fc)C+(fc) = -i"/k ■ uo{k). 
The solution C± (k) to the above system can be formally written as 

_ A+(fc)?io(fc) +i-fk-uo{k) _ -i-fk ■ uo{k) - X^{k)ho{k) 

^-^^>- A+(fc)-A_(fc) ' ^+W- A+(fc)-A_(fc) 

Therefore, n is solved as 

/O.N ~, ,N A+e^-*-A_e^+*. ,,, e^+* - e^-* , .,,, 

(3.4) n{t,k}^ no(fc) ijk-uo{k). 

A^ — A— A_|- — A— 

Denote k = k/\k\ when k ^ 0. From ik ■ E = —-fi, one has k ■ E = -TTjn. Thus, 



7"' ^-'"^ ''"■= "- "" ~ 7 Ifel' 

7lfc| -^^^^ i^uiiipuLiiig -p^,' 



multiplying p.4p by ^7|t and computing f t|t(— ^7^ • uo{k)) = (ik ■ uo{k), gives 



k.E= 1±— f^^ k ■ Eo + — ^ pk ■ uo{k). 

A+ — A— A-f — A— 

To solve fc • ii, from (|2.2p . by taking the time derivative of the second equation and 
then using the first and third equations to replace dtn, dtE, respectively, one has 

dtii + -I'^kk -u + liikx B + I3^u + n\kfdtu = 0. 
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Further taking the dot product with k gives 

dUk ■ it) + fi\k\^dt{k ■ u) + {j^\k\^ + (3^)k • u = 0. 
In the same way as for n, one can solve the above second order ODE as 



k ■ u{t, k) = -^— ^ k ■ u(0, k) 



A, - A_ 



-k-dtu{Q,k). 



A+ - A_ 
Compute from initial data 

k ■ dtu{0, k) ^ k ■ {—^ikho — /3Eq — /i|fc| uo) = ~«7|fc| no — (3k ■ Eo — /i|fc| k ■ uq. 
Therefore, 



(3.5) k-u{t,k) 



X+-\- 



-(-i7|fcpno) 



A+e^-* 



gA+t _ gA_t 



A+-A_ 



A+ - A_ A+ - A_ 

Pk-Eo). 



-(-^|fcp) 



k ■ iiQ 



Notice that from the indicate equation p.2p . 

\+ + \^ = -n\kf, X+\-=-f^\kf + p^, 
which imply 



A+e 



A_t 



"^-f 2n A+e^+*-A_e^-* 



-(-A#-n = 



A-f — A_ A-|- — A— A^ — A_ 

Applying the above identity into ()3.5p and dividing it by |A:|, one has 



k ■ u(t, k) 



A+-A_ 



(-i7|fc|7io) + 



A+e^+* - A- 
A+-A- 



,A_t 



-k ■ Uq 



gA+t _ gA_t 



-P~k -Eo). 



A+ - A_ 
In a word, one has proved 

Proposition 3.1. Let U = [n, u, E, B] satisfy the system p.l|) /or allt > 0, x <E M.^ 
with initial data U\t=o = Uq = [no,uo, Eq, Bq\. Denote k = k/\k\ when k ^0. Set 
u\\ by My = kk ■ u, and likewise for E\\ . Then, 

n \ I fiQ 

u\\ =A/' uo,|| 

except a finite number of values of \k\, where the matrix M ~ A/'(i, k)Yx7 is given 

by 



/A4 



,A_t 



Af 



-A_e^+ 



A+t 



A+-A- 



-ijk' 



A+-A_ 







\ 



gA+t_gA_t A-^e^+* - A_e-^-* 

-i^k- 



A+- 


-A_ 




A+- 


-A_ 


13 


-'A, 


-A_ 


-is 







/3 


gA+t 

A+ 


- A_ 


-t 

-h 


A+e^-* - 
A+- 


A_e^ 
A_ 


-I3 
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and X± = X±{k) given by (|3.3p denote the roots of the indicate equation (|3.2p . 

3.2. Green's function of the electromagnetic part. By applying the operator 
kx to equations of dtu, dtE, dtB in ()2.2p . one has 





dtk X u + pk X E + ^\k\'^k x u = 0, 
dtk X E — ik X (k x B) — j3k x u ~ Q, 
dtkx B + ik X {kx E)^0. 

Define 

Ml = Mi{t, k) = UA_{t, k) ~ —k X k X u, 
M2 = M2{t, k) = Ej_{t, k) = -kxkx E, 
M3 = Msit, k) = B_i_{t, k)^ -kxkxB. 

Then, Mi, M2, Afs satisfy the first-order ODE system 

"9tMi+M|A:pMi+/3M2 = 0, 
dtM2 -ikxMs- /3Mi = 0, 
^ dtM3 + ikx M2 = 0, 

namely, 

Ai|fc|2 ^ 

-/3 -ikx I ( A/2 I =0. 
ikx 

To simplify the computation of solving the above system, we take change of variables 

e = ^fc, r = f3t, 

and thus Mi, M2, M3 satisfy the rescaled ODE system 

/A/i\ /a|eP 1 

(3.6) dt\ Ahl + l -I -j^x 

VA/3/ V i^x 

with a = /i/3. 

Similar in [10], one can derive from p.6p a third order ODE satisfied by A/2- In 
fact, compute 9^ A/2 by replacing drMi, drM^ to get 

(3.7) dlM2 + (I^P + l)A/2 = -a|^pA/i. 

By taking one more r derivative of the above equation and then replacing d-rMi, 
it follows 

(3.8) dlKh + (leP + l)a.A/2 - a|epA/2 - a''\i\Hdi. 
Multiplying dXT]) by aj^p and adding it to ^^ yields 

dlM2 + a\£,\''dlM2 + (ICP + l)9rA/2 + alCI^A/2 = 0. 
Then, A/2 can be solved as 

A/2(t,0 = Cie^i- + C2e^^^ + C^e^'^ 
where \j = Aj(^), j ~ 1, 2, 3, denote the roots of the indicate equation 

(3.9) z^ + a\^\h^ + {\^\^ + l)z + a\^\*^0, zeC. 
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Here and hereafter we have excluded a finite number of values of |^| such that 
equation p.9|) has possible multiple real roots. Moreover, from equations of Mi 
and M3 in (|3.6|) . Mi(r, ^) and M3(t, ^) can be also solved as 

(3.10)Mi(r,^) = A/i(0,T)e-"l«l'''- / e-"l«l'(^--')M2(s,0rfs 

Jo 

3 



,=1 ^0 



e 



3X3 

and 

(3.11) M3(r,0 = M3(0,e) + / -«C X M2(s,0rfs 



= A/3(0,0-iCx / VCje^^'^ds 

3 3 

One can determine Cj = Cj{£,) in terms of Mj{0,^) as follows. The initial 
condition of p.6p implies 



I 1 i\ fcA ( M2(o,o 

Al A2A3 C2 = 9rA/2(0,C) 

.A?AlAi; \C^) \dlNh{Q,0, 

On the other hand, from equations of dr^h and d^A'l2 in p.6p and p.7p . 

M2(0,0 \ / 1 \ /Mi(0,O' 

9.M2(0,0 =1 z^x M2(0,e) 

_a^M2(o,o/ V-«i^i'-(i^i' + i) ; V^^^3(o,e). 

Combining the above two equations gives 

1 1 ly' f 1 \ /Mi(o,o 

Al A2A3 I I 1 i^x \ \ M2(0,0 

^1 A2 A3 
One can compute 

-1 




XfXlXV V-a|ep-(|?P + l) / VM3(0,0, 



111^ 1 

Al A2 A3 



^2 ^2 ^3 / (Al — A2)(A2 — A3)(A3 — Al) 



-A2A3(A2 — A3) (A2 + A3)(A2 — A3) — (A2 — A3) 

-A3Ai(A3-Ai) (A3 + Ai)(A3-Ai) -(A3 - Ai) 
^-AiA2(Ai-A2) (Ai+A2)(Ai-A2) -(Ai-A2)^ 
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Thus, 

(3.12) 

with C given by 




•Mi(0,C)~ 
= C| A/2(0,e) 
M3(0,C), 



(Ai-A2)(A2-A3)(A3-Ai) 



1, 



'(A2 + A3+a|en(A2-A3) (|eP + l-A2A3)(A2-A3) (A^ - Ai^x 

(A3 + Ai+a|en(A3-Ai) (|eP + l-A3Ai)(A3-Ai) {\l-\\)i^x 

_(Ai + A2 + a|^n(Ai-A2) (|eP + l-AiA2)(Ai-A2) (A? - A^^x 

Notice that 

3 

J|(z - Xj) = z^ - (Ai + A2 + A3)z2 + (A1A2 + A2A3 + A3Ai)2; - A1A2A3. 

Comparing the above with the indicate equation p.9p shows 

Ai + A2 + A3 = -a|^p, 

A1A2 + A2A3 + A3A1 = \i\^ 

A1A2A3 = -a|e|^ 
Then, one can rewrite C as 

^^■^^^ ^" (Ai-A2)(A2-A3)(A3-Ai) 

/-Ai(A2-A3) Ai(A2-A2) {\1-XI)i£,-k\ 

^A2(A3-Ai) A2(A2-A?) {\l-Xl)i^x 

V-A3(Ai-A2) A3(A?-Al) (A?-A2)z^x/ 

Basing on the formula p.l2p of the coefficients Ci(^) in terms of initial data A'/i(0, ^), 
it is a direct and long computation to verify the following lemma, and its proof is 
omitted for simplicity. 



Lemma 3.1. 



A'f3(0,O + *Cx^y^ = 0, 



J = l 



A/i(0,r) + ^ 



Cj 



j^,^j + 4i? 



From the above lemma, Mi(r, ^) and Mz{t,^) respectively given in p.lOp and 
p. lip are reduced to 

3 „ 3 



C, 



3 

3=1 ^ 



- = E 



c. 



^ 2-^v- 



+i 



X, 
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Thus, we have solved Mj{t, £_), j = 1, 2, 3, as 

■Afi(r,C)\ /Ci 

M2(r,C) Us C2 

,M3(r,C)/ V^3. 



with B defined by 

(3.14) B = 

Set 
Then, 



/ A2+A3 



pAlT 



A3 + A1 
0A2T 



e^3" \ 

>M+A2 \ 



oAsT 



-^^x^T -*^xV -*^xV/ 
7W =SC. 



/Mi(t,C)\ /Mi(0,^)\ 

M2(T,e) =>( M2(0,e) . 

VM3(r,e); VAf3(o,e)y 

In a summary, one has proved 

Proposition 3.2. Let U = [n, u, E, B] satisfy the system (|2.ip for allt > 0, x G R^ 
wit/i initial data U\t=o — Uq ~ [no, uo, Eq, Bo]. Denote k ~ k/\k\ when fc ^ 0. Set 

-uj_ = — fc X (fc X u) = (I3 — fc (g) k)u, 

and likewise for E± and B± . Then, 



u± 



k. 



uo.± 



Ei_ =M{pt,-) £;o,± 

except a finite number of values of \k\. Here, in terms of variables t = /3t and 
^ = ^k, the matrix M. = A^(t, ^)gx9 is equal to BC with B, C defined in (|3.14p and 
(|3.13p . respectively. Moreover, the elements of Ai{T,S,) take the form of 



^ ^^^^Mi)^f^ + e^--(-A2)^t^ + e^3r(_A3)A^ 
" (Ai-A2)(A2-A3)(A3-Ai) '' 

_ e^^-Ai(Ai - A^) + e^--A2(Ag - Af ) + e^^^X^{\\ - Aj) 
'' (Ai-A2)(A2-A3)(A3-Ai) "' 

'' (Ai-A2)(A2-A3)(A3-Ai) ' 



e^^"Ai(A2 - A3) + e^^-A2(A3 - Ai) + p>^-A3(Ai - A2)., 
-'^IS — — TTT"; TTT"; r^ -^3, 



Mi3 = 



(Ai-A2)(A2-A3)(A3-Ai) 

e^'^{X2 - A3) + e^^-(A3 - Ai) + e^^^^X^ - A2) 

(Ai-A2)(A2-A3)(A3-Ai) 



i^x, 



e^^^Xl Aj) + e^--(Ai - Xj) + e^^^Xj ~ Aj) , 
'' (Ai-A2)(A2-A3)(A3-Ai) ^ ' 
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M21 = -M12, M31 = Mi3, M32 = -M23, 
where Xj = Aj(^), j = 1,2,3, denote the roots of the indicate equation 

Now, by combining two representations of the Green's functions to the fluid 
part and the electromagnetic part obtained in Proposition 13. II and Proposition 13.21 
respectively, we have 

Theorem 3.1. Let U = [n,u,E,B] satisfy the system (|2.ip for all t > 0, x e M^ 
with initial data U\t=o = Uq = [nQ,UQ, Eq, Bq]. Assume that the solution U is 
written as 

U = G*Uo 
for the Green's function G = G{t, x). Then, under the decomposition 

with 

/ " \ f ^"\ / M 



c/(i 



Ul 

E\ 
\0 J 



, u^ 



(1) 



Wo, 
-^0, 

V J 



, (7(2) 



Uj_ 



, U^ 



2) 



/ \ 

Eo.± 



there are two Green's functions G^^' , G"^ defined by 



GW(t,fc)=(:^^i^ 







G(2)(i,fc) = 



-M(/3i,S) 



except a finite number of values of \k\, such that 

[/(i) === G^^) * C/($^'\ j = l,2. 

4. Solve a third order ODE 
In this section we study the property of roots of the indicate equation 

z' + a|^pz2 + (|^|2 + 1)^ + „|^|4 ^0, zee. 

Here, a > is a constant. In what follows, set 

r:=|eP>0, 

and 

g{z) := z + arz" + (r + l)z + ar . 

4.1. Elementary properties. We list some elementary properties for roots of 
g{z) = as follows. 

Property 1. Consider the existence of a real root for g{z) = with 2 e M. Check 
that 

• g{Q) = ar^ > 0; 

• g{—ar) = {—ar)^ + a^r^ — ar{r + 1) + ar^ = —ar < 0; 

• g'(z) = 3^2 + 2arz + (r + 1) > r + 1 for z > 0; 

• g'(z) = z^ + 2z(z + a?') + (r + 1) > r + 1 for z < -ar; 

• g{z) is strictly increasing over {z < —ar} U {z > 0}. 
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The above observation implies that equation g{z) = has at least one real root 
denoted by cr = a{r) G M, which satisfies —ar < a(r) < for all r > 0. Moreover, 
suppose that a = <j{r) e {—ar, 0) is a real root of g(z) = 0. Then 

g{(j) = cr^(o- + ar) + a{r + 1) + ar"^ = 0. 

Since cr^(cr + ar) > 0, then cr(r + 1) + ar^ < holds, i.e. 



cr < -- 



r + 1 
Therefore, any real root a — <T{r) of g{z) = actually satisfies 

„^2 



-ar < a{r) < 



r + 1' 



Vr >0. 



Property 2. Let r > i. Then —ar H — - G (— ar, 0). Compute 

g{-ar -\ ) = {-ar -\ )^ + ar{-ar H )^ + (r + l)(-ar H ) + ar'^ 



to obtain 



Notice 



1 N 1 
g(-arH ) = - 

ar a 



r ar ar 



0<^<l 

Then, whenever 
one has 



sup _ l + _ l-_ =_ 1 + —) 1- — 

1 ,1 ar ar ar 



VS V3 \/3 3\/3' 

3\/3 



3V^ 



a > 0, i.e. a < 



2 ' 



g{-ar-\ ) > 0. 

ar 

Thus, it follows that if a < ^^, then whenever r > i, equation g(z) = has at 
least one real root, denoted still by cr = cr(r), satisfying 

—ar < (7{r) < —ar -\ , Vr > — . 

ar a 

Property 3. Let a G {—ar, 0) be a real root of g{z) — 0. Consider the sign of 

g'{(7) = 3cr^ + 2ara + (1 + r) 
by the following two cases. 

Case i. Write 

g'(cr) = 3cr^ + r(2acr + !) + !. 
If2acr + 1>0, i.e. cr > -^, then ^'(ct) > 1. 
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Case ii. Using g{a) = 0, compute 

[3cr2+r(2acr + l) + l]o- 



g'i^) 



a 

.3 , o„^^2 



3(T'^ + 2ara^ + (1 + r)(T 



_ -3[ara^ + (1 + t')^' + ar^] + 2ara^ + (1 + r)a 

a 
-ara^ - 2(1 + r)G - iar^ 
a 
Write 

—ara" — 2(1 + r)(j — iar = —ar{a'' — 2) — 2(1 + r)(o' + ar) — ar . 

If cr2 _ 2 > 0, i.e. cr < -\/2, then g'(cr) > 0. 

Let —:j^< — \/2, i.e. a < ^. Then, from case i and case ii, one has g'{a) > 0. 
This implies that if a < -^ then equation g{z) = has one and only one real root 
a = (T{r). Moreover, the real root a{r) is strictly decreasing in ?- if a < -^. In fact. 



since 



then 



that is, 



From 



one can see 



g{<7) = <T^ + ara^ + {r + l)a + ar^ = 0, 

,g'((T)-^ + (acr^ + (y + ^ar) = 0, 
ar 

da —acr^ — a — 2ar 
dr 9'{'^) 

-aa^ ~ a — 2ar = —aa^ — {a + ar) — ar, —ar < a < 0, 



-aa — a — 2ar < 0. 



Therefore, 

da 

— < 0, Vr >0. 

dr 

Property 4- Let a = a{r) G (— ar, 0) be one real root of .g(z) = 0. Then 

g{z) = {z- a)[z'^ + {a + ar)z + (cr^ + ara + r + 1)] = 0. 

The rest two roots, denoted by x±i satisfy the equation 

z^ + {a + ar)z + (cr^ + ara + r + 1) = 0. 
Thus x± is solved as 

1. .1 



X± = --{a + ar) ± -z\/(sy + ar)^ — \(a'^ + ara + r + 1). 
To the end, we set (/J to be a function of a given by 

(^ = 3o-^ + 2ara - a^r"^ + 4(r + 1). 
It is straightforward to check 

1. ^ , 1 . 
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Claim: If a < ^ then (p > for all r > 0, and moreover 

-ir^+4r<Lp<4{l + r), Vr>-— . 
a r 2a 

Proof of Claim. Let a < ^. From Properties 1-3, a G {—ar, 0) is one and only one 
real root oi g{z) — and satisfies 

—ar < a(r) < —ar H , Vr > — . 

ar a 

Write 

(^ = 3cr^ + 2ar{a + ar) + r{4 - 3a^r) + 4. 

Then, whenever 4 — 3a^r > 0, i.e. r < ^Ar, o9 > 4 holds. On the other hand, 

2a ■ 



whenever r > - and — < ^W-, i-e. r > ^. one has 



> 3(7^ + 2ara - a^r^ > 3(-ar H )'^ + 2ar{-ar H ) - aV^ 

ar ar 

= 3{a^r^ - 2 + ^— ) - 2a2r3 + 2 - a^r^ 



which implies 

3 

4(1 + r) > 3cr^ + 2ara - a^r^ + 4(1 + r) = (^ > ^5-^ + 4r. 

Since ^ < gl^- is equivalent with a < ^^ and the inequality ^-^ > -^ holds, 
93 > holds true for any r > 0. Claim follows. 

Therefore, when a < ■^, equation g(z) = has one real root a and two non-real 
conjugate roots x± with 

4.2. Representation of roots. We shall consider a representation of all roots 
2; € C of g{z) ~ 0. For this moment, although we have known that there is at least 
one real root, it is not clear whether there exists an exact condition to assure its 
uniqueness. Recall the equation 

g{z) = z^ + C2z'^ + c\z + Co = 0, C2 = ar, ci = r + 1, cq = ar'^ . 

Define 

S ^2c\~ 9c2Ci + 27co = 2a^r^ - 9ar{r + 1) + 27ar'^, 
R^S^' - 4{cl - icif 

= [2a^r^ - 9ar(r + 1) + 27ar'^f - 4\a^r^ - 3(r + \)f . 

Then, there are two distinct cases: 

• i? > in which case there are one real root and two non-real conjugate 
roots; 

• i? < in which case there are three real roots. 
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The roots of g{z) = can be written as 

1 




^ = —— + 7:v3i, Lo = — - — —vSi. 



sJS + VR _ 3JS + i^f^R _ ^S^ -R^^^Q w _ , 2_o u/2„f 



= (^-^)^^''e* = (c^-3ci)^^^e* 



S - VR 3 S - i\/^R S'^ - R -^/Q ^i» / 2 o ^l/2 -^ 



V 2 V 2 '4 

where 

~ arccos ■ 



Notice C2 - 3ci > if i? < 0, that is, 

R^S^ - 4(4 - 3cif < 0. 
In this way, it is easy to see that for R < 0, 



3 S+y^. 3 S-yfR 

are a pair of conjugate complex numbers, and all roots A,;, i = 1, 2, 3, are real . 
Let's compute R. In fact, 

R = 27[4(1 + rf + aV2(8r2 - 20r - 1) + Aa^r% 

It is straightforward to observe that whenever r > is sufficiently small or large, 
R> Q holds true. In addition, since R can be also written as 

R = 27[4 + 12r + (12 - a^)?'^ + (4 - 20a^)r^ + 8a V^ + Aa^r^ 

then, as long as a < -j=, one also has i? > for all ?- > 0. 

4.3. Asymptotic of roots at and cx). Suppose that there are e — e{a) > 0, 
L = L{a) > such that R > over r < e(a) and r > L{a). Then, over {r < 
e(a)} U {r > L(a)}, equation g(z) = has only one real root a and two non-real 
conjugate complex roots x±j and thus we let 

Ai = (T, A2,3 = X±- 
Case when r < e{a): Xi = a = cr(r) satisfies 

g{a) = (7^ + ara + (1 + r)a + ar^ = 0. 
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Set 



a{r) = J2 



amr 

m=0 

as r — > 0. One can compute 

ao = fli = 0, 02 = ^o■^ 03 ~ o,^ 04 = — fl, 0,5 = a(l — a ). 
Thus, the real root <t can be written as 

cr(r) = -ar"^ + ar^ - or'* + a(l - a2)r^ H = -ar^(l - r + O(l)r^) 

as r — ?> 0. RecaU that two non-real roots x± take the form of 

(T + ar 1 
X±= ^±2^*- 

Then, it follows 

— ar^ + ar'^ — 0(l)r^ + ar ar(l — r + O(l)r^) 
«X+ == 7, = 7, > 

9x+ = -v/3a2 + 2ara - a'^r^ + 4(1 + r) 

= iv/3a2r4(l -r + 0{l)r^)^ - 2a'^r^{l -r + C'(l)r2) - a^r^ + 4(1 + r) 
= -^-a^r'^ - 0(l)r3 + 4(1 + r) 



= Vl + r-0(l)r2, 
which implies 

ar(l-r + 0(l)r2) , ^, , , 

X± = ^ 2 ^ ^^/\^r-0{\)T■^ 

as 1 — ^ 0. 

Case when r > L{a): Notice that 

g{(T) = (T^ + ara + (1 + r)(j + ar^ = 

is equivalent with 



,(T,o ,<T,, 1 + 7- cr a 



Set 



as r — >■ c». One can check 



oo 



r 



1 1 1,1, 1,1, 

ao = -a, ai =0, 02 = -, 03 = ^, 04 = ^(1 + -), 05 = =-(4 + -). 

a a"^ a'' a a° a 



Then, the real root a can be written as 

c l,lol,, l./l, 1 

- = -a+ -r-^ - —r-^ + —(1 + -)r-^ - —(1 + - 
r a o'' a'' a w^ a 



that is, 



1 /. 1 0(1), 

a = -ar -\ (1 5- + — V^) 

ar a'^r r^ 



as r — ^ 00. For non-real roots, let's further compute 

(T + ar I , 1 0(1) 

— (1 1 

2ar a^r r 



^x± = -^^ = -:^(i-::^ + ^), 
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and 



(^ = 3(7^ + 2ar(j - a^r^ + 4(1 + r) 

~ 3cr(cr + ar) — ar{(7 + ar) + 4(1 + r) 

= r( a){(j + ar) + 4(1 + r) 

r 

3 ,, 1 0(1), 



-4 
-4(1- 






1 



_^ , o(i) 
0(1) 



1,1 -^ + 20),+ 4(1 



,,2 



(l-i + ^)+4(l + -) 



a^r 






1 0(1 



2 +^)' + 4(l + -) 



= 4r + 0(l)r" 
Therefore, one has 

X± 
as r — > oo. 



1 — ^^j, j,^ 



1 

2^r" 



We conclude this section with a remark. The asymptotic behavior of roots A^ (r) 
at origin and infinity exactly determines the same thing for the Fourier transform M 
of the Green's function of the electromagnetic part. However, M. is not defined at a 
finite number of values of r = |^p = I^P//3 where i? = 0, that is, possible multiple 
real roots occur. Therefore, to estimate M. over the whole frequency domain, one 
has to characterize the behavior of all roots \i(r) in the domain {e(a) < r < i(a)}, 
particularly over some small neighborhood of a finite number of points such that 
i? = 0; we shall clarify this later on. Notice that this issue is different from that 
in [To] where the situation is simpler due to the fact that all roots are far from the 
imaginary axis over {e{a) < r < L{a)} and multiple roots do not happen. 

5. POINTWISE ESTIMATE ON THE GREEN'S FUNCTION 

In this section, we consider the pointwise estimates on Af and Ai which arc the 
Fourier transforms of the Green's function to the fluid part and electromagnetic 
part, respectively. The most difficult part is the analysis of A^ and M in the finite 
domain far from origin and infinity, which turn out to decay in time exponentially. 
We then use these estimates to derive the corresponding pointwise estimates on the 
Fourier transform of solutions to the linearized Navier-Stokes-Maxwell system (|2.ip , 
and further obtain the time decay property of each component in the solution. 

5.1. Pointwise bound on the fluid part. Recall that X± satisfies the indicate 
equation 

z2 + M|fcpz + (7'|fc|'+/3')=0, 
and 

A± = -^^l\k\' ± iy^, i; := ^i'\k\^ - 4(72|fc|2 +/32). 

Divide the whole space M'^ by three subdomains as 

R^ =©oUDiUDoo, 
IDo = {|fc|<e}, ©i={e< |fc| <i}, Doc = {|fc|>L}, 
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and further set 

©r = {e < \k\ <L,ip< 0}, n^ = {e< \k\ < L, V > 0}. 

Case when \k\ is near 0; In this case, as |fc| — > 0, 

A± = -^M|A;|2±(/3 + 0(l)|fc|2)^ 
and 

1 2 \ \ ^,2 1 7^ 1 2 1 /j2 



Notice that in Afii,Af33, 

A+e^-*-A_e^+* _ A+A_ e^+* e^-* 7^|fcp + ;3%^e^+* 



A_|- — A_ A+ — A_ A+ A_ -^Z— -0 A-f 

Then, A/iijA/as can be estimated as 

|Mi| + |AA33|<e-5M|fe|^t. 

In the similar way, it holds that 

|A/'22| + |A/'23| < e-^^l*'-!'*, IM2I < |fc|e-^A'l'=l'*. 

Case u;/ien |A:| is near 00; In this case, one can compute that as |fc| — >■ 00, 



= - 


4a#P 


±^M|A:P(1 


/32 ' 


k\-'-Oil)\k\ 


-% 




which implies 














A_ = - 


-l^^\k\' 


1 ,,,2 


/32 + 


o{i)\k\-' = - 


-0(l)|fcp, 


A+ = - 


-^A#-P 


+ ^M|A:P- 


A«7 

/32 


o{i)\k\-'^- 


A«7 

/32 


-0(1)1^1-2 


Therefore, one 


has 












Wu 


1 + IAA33 


l< 1^+^- 
1- IA+-A 


-1 1^^+ 












/7'|fcp+/3Y' 


-{^+o(i)\k\-^) 


It 

e 


-0(l)|fe| = t 
-i/'iM;^I2 '' 



Vi^ ' i|^ + o(i)|fc|-2 o(i)|/cp 



< p-0(l)t I J_„~0{l)\k\h 

^ \k\^ 
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and 

< 1 ((^ + 0(l)|fcr2)^-(-^+0(i)|.|-)* ^ o(i)|fc|2^-0(i)|.|^*) 

In the similar way, one can also obtain 

|M2|< JL_(eA+* + e^-*) 

< J_e-o(i)t + J_e-o(i)|fept 



|fc| 



and 



|A/--l^^e-«^^^*+ul 



|fcP 



-o(i)|fc|^t 



Case w/ien |fc| is Far from and oo; In this case, A+ tends to A_ as |fc| is properly 
chosen such that 0^0. Still consider first the estimate on A/ii, A/ss. Notice 

A+e^-* - A_e^+* . . [^ d -^^ 



= -A+A_ / — ( — )|A=(l-s)A_+sA+rfs, 



A+ - A- Jo 

where 

_d_ e^ {t\ - l)f/* 

'dX^~'^ A2 ■ 
Then, it follows 

lAAnl + IAA33I < |A^A_| /' ^(i_i)^_±i^e*^«-^)-+-^>ds. 

We make further estimates over two subdomains D]^ and D^ as follows. 
Over D^ = {e. < \k\ < L ^ ij; > {f\: A± < are real, which shows 

(l-s)A_ +.sA+ < 0, 0< s < 1. 
Thus, one has 

lAAnl + IA/-33I < /' '"^'"'fn" 'I'- +/^+' + ^ e*-^""^{(^-)>-^-^^>d. 
Jo mm|(l-.s)A_ +sA+|^ 

<(l + i)e-o(i)*<e-oW* 
where max, min are taken over {fc G Dj*", < s < !}• 
Over Dj^ = {e < |fc| < -^,"0 < 0}: A± arc complex conjugate, which shows 

^{{\ - s)A_ + s\+] = SRA± = -i/ilfcP < 0. 
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Then, in the similar way, one has 



IMil + IA/-33I < /' ^=|^7-%i^^^e*^«^-)^-+^^+>d. 



jQ min|3f{{(l-s)A_+sA+}|2 

<(l + t)e-0(i)*<e-0W*. 

where max, min are taken over {fee Dj^, < s < 1}. 

Therefore, from both cases above, over Di = D^ U D^, it holds 

|AAn| + |AA33|<e"°W*. 

In the completely same way, over I?i, 

|AA22| + |M2| + |AA23|<e-o(i)*. 

In sum, we have proved 

Proposition 5.1. Let the matrix M = Af{t,k)jx7 be defined in Proposition \3.1 
The elements of M have the following upper hounds for pointwise t > and k E M.^ 
Over Do, 

( 1 \k\d^ \ 
\N\ < |fc|d3 I3 I3 e-^^l'^-l'*. 

V I3 hj 

Over Doo, 

1 1^1-14 \ 

lA^I < I \k\-'d, |fc|-2l3 |fc|-2l3 e-O(i)* 

|fc|-% I3 / 

|-2 lfcl-14 \ 

\k\-^d, I3 1^1-% e-Od)!^!^*. 



\k\-n^ |fc|-2l3 



OverDi, 



1 4 

\Af\ < I d, I3 I3 I e-«(i)* 

I3 I3, 

Here, for simplicity, \J\f\ means the matrix corresponding to J\f with each element 
taken the absolute value, d^ denotes the column vector (1,1,1)"^. 
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5.2. Pointwise bound on the electromagnetic part. In order to estimate the 
pointwise bound on A4, let's rewrite its elements as 



^ Ai 

' (Ai-A2)(Ai-A3)(A2 + A3)' 
1 
(Ai+A2)(Ai + A3)(A2-A3) 

Ai(A2+A3) 
(Ai-A2)(Ai-A3) 
(Ai+A2)(Ai + A3) r A 



Ml 



(Ai-A2)(Ai-A3) 



.Air 



A2(Ai + A2)^A2r _ A3(Ai + X3) ^\:,r 



Ai — A2 



Ai — A3 



A2 — A3 



2_^X.r 



XI -xr A? -a: 



^e^3. 



^ ^ ?-(A2+A3) Air 

"'•' Ai(Ai-A2)(Ai-A3) 
r(Ai + A2)(Ai + A3) 



At — A3 



1 



„^2T 



0A3T 



A2(A?-Ai) A3(Af-A2) 



and 



M12 = 



Ai 



(Ai-A2)(Ai-A3) 



.Air 



1 



A2 — A3 
Mi3 = {- 



X 



2 gA.r _ 



Ai — A2 
1 



^^gA3r 



Ai — A3 



(Ai-A2)(Ai-A3) 



=Air 



1 



A9 — A3 



1 



Ai — A2 



oA2r 



.Asr 



Ai — A3 



i^x, 



M23 = i - 



A2 + A3 



.Air 



(Ai-A2)(Ai-A3) 
(Ai + A2)(Ai + A3) 
A2 — A3 



.A2r 



.Asr 



A-^ A9 A-^ A3 



»Cx, 



and 



M21 = -M12, M31 = Mi3, M32 = -M23- 

Similarly as for considering A/", let us divide the whole space R"^ by three subdomains 
as 

M^ =DoUDiUDoo, 
Do = {r < e}, Di = {e < r < L}, Doo = {r > L}. 

Here, we still used Do,ID'i,]D'oo for simplicity of notations without any confusion. 

Case when r is near 0, i.e. r G Do- In this case, Ai, i = 1, 2, 3, have the following 
asymptotic behavior 

Ai = cr(r) = -ar^{l - r + 0{l)r^), 

ar — ar'^ + 0(l)r^ , , , ^ 



A3 =X- =X+, 
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as r tends to 0. This also implies 



Al — A2,3 „ 



\ar - 0{iy T Wl + r-0{l)r-^, |Ai - \2,z\ = 0(1), 



1 



Al + A2,3 = "2"'' " ^(^)''' ^ ^Vl + ^-0(l)r^ |Ai + A2,3| - 0(1), 
A2 + A3 = -{ar - ar^ + 0{iy), IA2 + A3I = 0(l)r, 
A2- A3 = 21^/1 + ^-0(l)r2, |A2-A3|=0(l)ri/2. 

By using the above computations, it is straightforward to make estimates on M as 
r tends to 0: 

|-M33|<e-°('''^'" + re-o<'^'■^ 
\M23\ < r3/2e-o(i)'-'- + ri/2e-o(i)rr^ 



Case when r is near cx3, i.e. r G Doo- In this case, one has 

1 /, 1 0(1), 
Al = (T = -ar + — 1 - ^ + ^-Z ^ 

2ar a^r r 

A3 = X- = X+ 



and 



^2 = X+ = - — {l- — + ^) + i^r + 0{l)r'\ 



as r tends to 00. This further implies 

0(1 



Al - A2 = -ar + ^ - i^r + 0{l)r-\ |Ai - A2I = 0(l)r 
r 

o(i) 



Al + A2 = -ar + —^ + iyjr + 0(l)r-i, |Ai + A2H 0(l)r, 
r 

A2-A3 = 2zVr + 0(l)r-i, IA2 - A3I = 0(l)ri/2, 

A2 + A3 = --(1 - ^ + ^), IA2 + A3I = 0(\)r-\ 
ar a r r^ 

Then, as in the case when r is near 0, it is also straightforward to obtain the 
estimates on A^ 



i-M33i< V°^^^'^^ + e-'^^ 
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and 



|Ali2|<-e-oW- + ie- 



1 o(i) , 
— ( 
r r 






as r tends to 



Case when r > is far from and oo, i.e. r € Di. In this case, we shall consider 
the estimates on Ai by dividing Bi = {e < r < L} into three parts. For that, 
denote A^^l as a set of all zeros of R over [e,L], that is, 

^e,L = {?^», l<i<io}-={e<r<L,R = 0}. 

Notice that from Property 1 in Section 14711 g{z) = can not have one real root 
with three multiples. Therefore, for any ri G A^^l, by the representation of roots in 
Section|421 g{z) = at r = r^ has one 1-multiple real root and the other different 2- 
multiple real root. Set a to be the 1-multiple real root, and x+ = X- the 2-multiple 
real roots, of g{z) = at r^. Then, it follows 

g'{a)\r=r, > 0, ip\r=ri = 0, 

where ip depending on a and r is given by 

(^ = 3(7^ + 2ara - a^r'^ + 4(r + 1). 
Define 

(Jo = o IpiP 9'{<j)\r=r, > 0. 

Z 1<2<20 

By continuity, for any r) > small enough, there are ei-ncighborhoods B{ri,€i) C 
[e, L] which do not intersect, such that 

min inf g' {<y) > (5o, max sup |(p| < rj. 

l<i<io B(ri,ei) l<i<io B{ri,ti) 

We now divide two cases to decide how to choose roots Ai, A2, A3 and further make 
some essential estimates. 

Case 1. For each 1 < « < ?0i over -B(ri, e^), we choose a depending on ri as before, 
and set 

Ai = cr, 

G + ar I 

A2,3 = X± = ^ — ± 2^ -'P- 



One can check 



which implies 



3(7 + ar 1 

Ai - A2,3 = ^ T 2 V ~'^' 



(Ai - A2)(Ai - A3) = i[(3(7 + arf - {-^)] = g' {a) > <5o > 0. 
Moreover, one can choose ry > small enough such that 



^»-|-\/l>«- 
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Then, it holds that 






and 



|3cr + ar\ 






Case 2. Over [e, LWul'L^ B(ri, e^), we choose a as the minimum (or maximum) real 
root of g{z) = 0. Then, by this choice, one has 

g'(a)>0, 1^1 >0, yre[e,R]\UtiB{n,e,). 

This implies 

(Ai-A2)(Ai-A3) = i[(3a + a7f-(-(^)]=g'(a)> min g'ia}>0, 



le,R]\uZ,B{n,e,) 



and 



min \ip\ > 0. 

[e,R]\uZ,B{r,,e,) 



We now decompose 



Di,- = {e < r < L, r ^ UtiB{r,,e,), R < 0}, 
Bi,+ = {e<r<L,r<^ \JtiB{n,€,), R> 0}, 
^1,0 = {e <r < L,r eUtiB{r,,€,)}. 
Notice that all domains Di,± and Bio are bounded and closed, and 

-R|di.+ > 0, i?|Di,_ < 0. 

In what follows we shall estimate Alu only; other elements of Al can be estimated 
in the completely same way. 



Over Di^q: Rewrite Mn as 










■^"^ (Ai-A2)(Ai-A3)(A2 + A3)'' ^' 




1 r^ d 


rA(Ai + A) ,1 
Ai-A 




{Xi + X2){Xi + Xs) Jo dX 


A=(l-e)A2+6/A3 


where 

d 

dX 


rA(Ai+A) ,1 
Ai-A ' 


tA(Ai + A) + Ai + (Ai + 2)A + A^ ^^ 
(Ai-A)2 " ■ 


Then, one has 


(5.1) \Mn 


'- |(Ai-A2)(Ai-A3)|-|(A2- 
1 


fA3)l 


|(Ai+A2)(Ai+A3)| 


Jo 


r|A(Ai+A)| + 
1- 


|Ai + (Ai+2)A 
\i-A|2 


+^\'«] 


A=( 


de. 

i-9)X2+ex3 



ddh 
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Notice 



Also, notice 



which implies 



A2 + A3 = -{a + ar) < 0, 
3?Ai = cr < 0, 
3fi{(l-^)A2 + 0A3}<O. 



a — ar I 

Ai + A2,3 = — ^ — ± 2 V ~'^' 



(A, + A2)(A, + A3)=(^-"'^)' + ^ 



4 
{a - arY + Scr^ + 2ar(T - a^r"^ + 4(r + 1) 



cr^ +r + 1 > 0. 



Consider the low bound of |Ai - A| at A == (1 - 6')A2 + 6* A3 for all < 6* < 1. 
Case when A2,3 are real: In this case, 

Ai - [(1 - e)x, + ex,] = ^^ - {20 1)^ 

By the definition of Di.o, 



|A,-[(1-^)A2 + ^A3]|>^^-^M> /^_ /|>0. 



Case when A2,3 are non-real and thus complex conjugate: In this case, 

Ai - [(1 - 6i)A2 + 61 A3] = Ai - 5RA2 + (261 - l)i3A2, 
which also implies 

|Ai - [(1 - 9)X2 + ex^w > |Ai - 5RA2I = ^^±^ - Vf ^ °- 

Therefore, whether A2,3 are real or not, 

min |Ai-[(l-6i)A2 + 6iA3]| >0. 

reDi.o,o<e<i 

Then, putting all estimates above into (|5.ip gives 

l-^iil < (1 + T)e-O(i)^ < e-o(i'^ r e ©1,0. 

Over Di.-(_; By noticing 

A2 + A3 = 23fiA2 = -(a + ar) < 0, A2 - A3 = 243A2 = iy^, ifi > 0, 
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it is straightforward to make the estimate 



2 |A2(Ai + A2)|^,5,,, 



< 



|(Ai + A2)(Ai + A3)|-V^ IA1-A2I 
max I All 
min |(Ai — A2)(Ai — As)] • min \a + ar\ 
2 



min|(Ai + A2)(Ai + As)] • \/miii Lp 
max|A2(Ai + A2)(Ai - As)] rmin -(■^+°'-) 
min|(Ai - A2)(Ai - As)] 

< p-0(l)r 

where max, min are taken over Di,+. 
Over Di,_; Since 

A2 + A3 = 23fiA2 = -(o- + ar) < 0, A2 - A3 = ^/^, 1^ < 0, 
it foUows in the same way that 

Thus, for 7M11, one has 

l-^iil ^ e-°(i)^, r e Pi = Di,o U X>i,+ U Vi,^. 

The above estimate also holds true for all other elements of M and details are 
omitted for simplicity. 

Then, in a summary, we have proved 



Proposition 5.2. Let the matrix AA = A4{t,k)Qx9 be defined in Provosition \3 .^ 
The elements of Ai have the following upper bounds for pointwise t > and /c G M'^. 
Over Bo, 

( r r2 ri/2\ / l 1 ^1/2 \ 

\M\< 7-2 r3 r3/2 e-°(i)'-'"+ 1 1 r^'^ e~^'^^>\ 

\^j,l/2 ^3/2 I J \^^l/2 ^1/2 ^ y 

Over Doo, 

|X| < r-i r-2 r-5/2 g-^^i)''^ + r"! 1 1 e"^^. 
\r-3/2 r-^/2 r--^ / \r-^ 1 I J 

OverDi, 

/111\ 
|7W| < 11 1 e-o(l)^ 

By combining Proposition 15.11 and Proposition 15.21 together with Theorem 13. 1[ 
we have the following 
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Theorem 5.1. Let U = [n,u,E,B] satisfy the system p.ip for all t > 0, x iz 
with initial data U\t=o ~ Uo ~ [no,uo, Eq, Bq]. Then, it holds that 

f\n\\ f\no\\ 



\u\ 

E 

\\B\J 



< (JuPP 



for all fc e R'^, where over Do, 



G"PP = 



/ 1 |fc| \ 

|fc| 1 1 |fc| 

1 1 |fc| 

\Q |A:||fc||A:|V 



-0{l)\k\H 



l"ol 




l^ol 




\\Bo\) 


/O 

0|fc| 

0|fc| 

\0 \k 


\ 

4 |fc|6 |fc|3 

\k? 1 1 



-o{i)\k\H. 



over . 



^upp _ 



/ 1 |fc|-i 0\ 

|fc|-l |fc|-2 |fc|-2 _ 

|fc|-2 1 "" 
\ 0/ 

/|fc|-2|fc|-i \ 

|fc|-l 1 |fc|-2 |fc|-3 
|fc|-2 |fc|-2 |fc|-5 
\ |fc|-3 |fc|-5 |fc|-6y 



o(i)t 



-0(l)|fc| = t 



/O \ 

0|fc|-4 |fc|-2 |fc|-2 

0|fc|-2 1 1 
Vo|fc|-2 1 1 / 



0(l)t 



,-o(i)t 



/I 1 00\ 

/^upp i i 1 1 

111 

Voiii/ 

5.3. Refined time-decay property. From Theorem 15.11 for any integer j > 0, 
one has the following estimates for t > 0, fc G K"^. For n, 

|.F(V%)| < (IfcHnol + |fcP+^|^o|)e-^l'=l'*lBo 
+ \k\mho\ + \uo\)e--'lB, 



-{\k\^\m\ + \kr'\uo\)e 






Here, for later use, we use hq = —-gik ■ Eq to replace no in the first term on the 
right-hand side so that 

|^(V%)| < i\k\^+'\Eo\ + |fcr^|iio|)e-^l'^-|'*lDo 
+ \k\\\m\ + \uo\)e~''lB, 
+ i\k\^\no\ + \k\^-^\uo\)e-''lo^. 
For u, 

\T{\/^u)\ < i\k\^+'\ho\ + Ik^uol + Ik^Eol + |fcp+i|Bo|)e-^l'=l'*loo 
+(|fcP+2|uo| + Ifc^iol + |fcP+i|So|)e-^l^l'*lDo 
+{\kr'\no\ + Ikmiiol + \kr^\Eo\ + |fcp-3|^o|)e-^*lD,uD^ 
+(|fcp-4|tio| + \kr^\Eo\ + \kr^\Bo\)e-i^lo^. 



and 



GREEN'S FUNCTION OF THE NAVIER-STOKES-MAXWELL SYSTEM 41 

For E a.nd B, 

l-^(V^'^)| < (lA^n^ol + l^n^ol + |fcP+^|So|)e-=l'=l'*loo 

+ (|fcP+4|zlo| + lA^r^l^ol + |fcP+'|So|)e-^l^-|'*lD„ 
+ {\k\^-^\uo\ + Ik^Eol + \kr'\Bo\)e-''la,uB^ 
+i\kr^\M + Ik^Eol + l^nSoDe" w^Id^, 

\n'^'B)\ < (|fcp"+i|^o| + \k\'+'\Eo\ + |fcP+2|So|)e-=l'=l''lDo 
+ (|fcP+i|zlo| + \k\'+'\Eo\ + |fc|-'"|So|)e-^l'=l'*lDo 
+ i\kr^\uo\ + \kr^\Eo\ + \kr'\Bo\)e-^*ln,uB^ 

+ {\kr^\uo\ + \k\^E„\ + \kmBo\)e~i^la^. 

Now, for brevity of presentation, we introduce a notation related to the rate index 
over the low-frequency domain as follows. For integer j > and 1 < p < q < oo, 
define 

Sij,P,q) = -^ + T:i 

2 2 p q 

which corresponds to the usual value in the case of the heat kernel in M.^. Then, by 
applying Lemma 12.11 to those estimates above, it is straightforward to obtain 

Theorem 5.2. Let U = [n,u,E,B] satisfy the system (|2.ip for allt> 0, x e M^ 
with initial data U\t=o = C/q = [?^0j Wq, i?o, Sq]. Let l<p,s<2<q<oo 
and let £ > 0. Then, for any integer j > 0, the following time decay property of 
U = [n, u, E, B] holds true: 

\\y'n\\L, < (l +t)-'^J+'-P^'^^\\Eo\\Lp + {I + t)-'<^^+'^P^^^\\uo\\Lp 
II V^wIIl, < (1 + f)-^(^+i'P'9) II [no, Bo]\\lp + (1 + ty'^^'P^"^ II K, Eo]\\l. 

+ (l + t)-^*(^+2'P-9)||uo||LP + (l + i)^^*'''+^'^''^ll^o|lLP 

+(l + t)-^*(^+i'P-9'||Bo||Lp 

+g-ct(j|ymb-l,.,9)„^j|^_, + llV^f^'-^'^^UolU^ 

+||V™(^-2'^'«)^o||l= + ||V™(^-='''''«)Bo||lO 

+il + t)-^\\V"'^i-*+'-'-'J^uo\\L~^ + \\V'^^'-^+'-''''^[Eo,Ba]\\L^), 

\\^'E\\L.<(l + t)-'^^^P''J^\\[uo,Eo]\\L. + {l + tr'^'+''''''^\\Bo\\L. 

+ (l + t)"^*''^'+^'^'''^||Bo||LP 

+e-ct(||vm0--2,.,g)^^||^^ + ||V™(^'^-'')-Bo||l^^ + ||V™(^"5-^'«)Bo||lO 
+ (l + f)-5(j|v™(-'"-2+^'-9)y^j|^^, + ||V™0+^^«>9)[i?o,Bo]||i.,) 
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+ (1 + f)-^*-W+i.P.9) WuoWl. + (1 + t)-^*(^+3^J'^9) i^EoWl. 
+ (l + f)-3^(^-P^?)||Bo||Lp 

/or any i > 0. 

For later use, we need the following result which is an immediate corollary from 
Theorem [ 



Corollary 5.1. Under the assumption of Theorem 1 5. 2\. one has 

ll"-|| < (l + i)^^l|["-0,'"o]||LinL2, 

11^*11 < (l + t)-^\no\\L^nL^ + Ci- + ty-^\\[uo, Eo]\\L^nL^ + {^ + tyhlBoh^nL^, 
\\E\\ < {1 + t)-'^\uo\\LinL- + (l + t)-'H\\En\\LinL- + \\y^En\\) 



+ il + t)-m\Ba\\L^nL-^ + \\V^Bo 



and 



\\B\\ < {l+t)--s\\uo\\LinL-^ + (1 + t)-i{\\Eo\\L^nL^ + \\'^''Eo\\) 

+ (l + f)-i(||So||LinL2 + ||VBo||: 
for any t > 0. Moreover, one also has 

\\VB\\<{l + t)-^\uo\\L^nL^- + il + t)-l{\\[Eo,Bo]\\L^nL^- + \\V''[Eo,Bo]\\) 
for any t > 0. 



6. Nonlinear asymptotic stability 

In this section we are prepared to prove Theorem 11.11 concerning the global 
existence, uniqueness and time decay rates of solutions to the reformulated Cauchy 
problem (|1.3p . (|1.4p . (|1.5p for the nonlinear Navier-Stokes-Maxwell system. We 
first prove the global existence and uniqueness by only deriving the key uniform- 
in-time a priori energy estimates, and then apply the time decay property of the 
linearized system together with the further high-order energy estimates as well as 
the bootstrap argument to obtain the desired time decay rates of solutions. 

6.1. Global existence. First of all, still using the notation [n, u, E, B] for simplic- 
ity, let us write the reformulated nonlinear system (|1.3p . (|1.4p as 

' dtn + 7V ■ u = hi, 
dtu + jWn + /3E - ^Au = /12 
dtE-V X B- Pu = ha, 
dtB + V X E^O, 

S/-E^--n, V-B = 0, 

7 



(6.1) 
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where V • /13 = — -/ii and the nonUnear source terms hi,h2, ft.3 behave as 

hi r^ \7 ■ {nu), 

h2 '^ u ■ Vu + riS/n + u x _B + nAu, 
/13 ^ nu. 
Initial data of the system is given by 

(6.2) n\t=o = no, u\t=o = uq, E\t=o = Eq, B\t=Q = Bq 
satisfying 

(6.3) W ■Ea = --no, V • Bo = 0. 

7 

Then, we have 

Theorem 6.1. Let N > A. Assume (|6.3p for initial data Uq := [jiqj wo,i*^0: Sq]. // 
||t/o!lH« is sufficiently small, then the Cauchy problem (|6.ip . (j6.2|l admits a unique 
global solution U = [n, u, i?, B] with 

C/GC([0,oo);i/^(R3)), 

n e L^iiO, 00); i/^(M3)), Vm e L^{{0, 00); ^^(M^)), 
V-B e L2((0,oo);iJ^-2(M3))^ y^B e L2((0,oo); iJ^-^^M^^), 
and 

(6.4) \\U(t)\\l. + f (||n(s)||^„ + \\\7u{s)\\l. 

+ !|VS(s)!Ih,v-2 + \\V^B{s)\\]j^^,)ds < C\\Uo\\hn 
for any t > 0. 

Proof. We only consider the proof of the following uniform-in-time a priori estimate. 
Define 

£N{U{t)) = \\[n,u,E,B]\\l^+Ki Y. (9"(7Vn).a"«) 

|q|<7V-1 

+Ki Y^ (a"V X E ■ d°V X u) 

\a\<N-2 

+K1K2 XI (a"(-VxB)-9"£:), 

l<|a|<Ar-2 

and 

2?A.(C/(t)) - |ln|||,« + !|Vt*!l^« + |lVii;f^„_. + ||V2i?||^„_3. 

Here, constants < ki,K2 ^ 1 are properly chosen as in (|2.12p . Then, for any 
smooth solution U to the Cauchy problem dtH), ([O]) over < i < T with T > 0, 
one has 

(6.5) ^£N{U{t)) + cDNiUit)) < C{\\n\\L^)8N{U{t)y/^VN{U{t)) 
at 



for any < t < 2^. As long as the above estimate is proved. Theorem 16.11 follows in 
the standard way by combining it with the local-in-time existence and uniqueness 
as well as the continuity argument. Therefore, in what follows we prove (j6.5p only, 
and other details are omitted for simplicity. 
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In fact, zero-order energy estimate implies 



= J2 i{d"'hi,d"n) + (9"/i2 • 6>°u) + {d^hs ■ d^E)). 

\a\<N 

In the same way as in the proof of Theorem I2.1[ from the first two equations of 
(|SJ|) and V ■ E = -^n, one has 

J- Yl (5"(7Vn),a"w)+c Y, / (7'|5"Vnp+/32|9"nndx 

|Q|<Af-l |a|<7V-l"^*^ 

\a\<N-l |Q|<Af-l 

By also using the curl of the second equation of (|6.ip . i.e. 

dtV X u + 13V X E~ ^AV x u = V x /i2 
together with the third equation of (|6.ip . it follows 



^ Y {d"V xE,d°'\7 xu)+c J2 [ \d°'VxE\'^dx 

|q|<JV-2 |Q|<Af-2 '^'^ 

<e y / \d"yxB\^dx+-\\Vu\\l, 



l<|Q|<Af-2' 

+ J2 (^"^ ^ /l3, 9"V XU)+ Y (^"^ ^ -^' ^"^ >< ^2> 

|a|<JV-2 |q|<7V-2 

for an arbitrary constant < e < 1- Moreover, the combination of the third and 
fourth equations in (|6.ip gives 

4- Y {d°'{-V xB),d°'E)+c y / |9"VxBPda; 
at ^-^ ^—' Jths 

y /" |9"V x£;|2dx + C y / |9" 



l<|Q|<Ar-2 l<|a|<JV-2 

< ^ / |9"V x£;|"dx + C y /" Wu?-dx 

1<|q|<JV-2^'^ l<|a|<Af-2' 

y {d-^V xB,d"h:i). 



l<|a|<JV-2 
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Then, by putting the above estimates together for properly chosen constants < 
Ki, K2, e ^ 1, one has 

(6.6) j^£NiU{t))+cVN{U{t)) 

< J2 {{d"hi,d°'n) + {d°'h2,d"u) + {d"h3,d°'E)) 

\a\<N 

+Ki Yl (a"(7V/ii),a"u) + Ki Yl (5"(7Vn),a";i2) 

|Q|<Ar-l |a|<Af-l 

\a\<N-2 |a|<JV-2 

l<\a\<N-2 

The rest is a long and standard procedure to verify by using the Leibniz formula 
and the Sobolev inequalities [TJ [55] that the whole right-hand term of the above 
inequality is bounded by C(||n||L~)fAr(C/(i))^/^I'Ar([/(i)), where C depends only 
on llnllioo; for simplicity, we also omit the proof's details of the rest part, cf. [TT| 
and [TD]. The proof of Theorem 16. II is complete. D 

6.2. Large-time behavior. From now, we suppose that all conditions in Theorem 
l6.1l hold and U = [n, u, E, B] is the obtained solution to the Cauchy problem (|6.ip . 
(|6.2p with the condition (|6.3p . In this subsection we devote ourselves to proving 
the time decay rate of the full energy ||C/(t)|l^iv or cquivalcntly £N{U{t)). For that 
purpose, define 

X{t)^ sup (l + s)3£Ar([/(s)), t>0. 
0<s<t 

In fact, we have 

Lemma 6.1. // j|C/o||Linff™+i is sufficiently small, then 

inff«+i. 



(6.7) supX(t) < C\\Uo\\l^nm 



t>o 
Proof. Under smallness assumption of ||C/o||if«! (|6.5I) implies 

(6.8) ^^EN{U{t)) + cPjv(t/(t)) < 

for any i > 0. This is the starting point to deduce ()6.7p . In fact, fix a constant 
e > small enough. Then, the further time weighted estimate on (j6.8p gives 

(1 + t)-^+'EN{JJ{t)) + c f {1 + s)i+'VN{U{s))ds 
Jo 

< £n{Uo) + (7 + e) / (1 + sri+'8N{U{s))ds. 
^ Jo 

Noticing 

£n{U) ^ |1(7|1^„ < Vn+i{U) + \\[u,E,B]\\' + \\VBf 

<CVN+iiU) + \\[u,E,B]f 
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and 

VN+i{U{s))ds<C£N+i{Uo), 
/o 

it follows that 

(6.9) (1 + t)i+'£N{U{t)) + / (1 + s)i+'VNiU{s))ds 





<C£n+i{Uo)+C f {l + s)-i+'\\[u,E,B]fds. 
Jo 



Due to Corollary 15. 11 one has 

\\B\\ < Cil +t)-i\\[uo,Eo,Bo\\L^nHS 

+C f (l+t-s)-t||/i2(s)||iinL2ds 



+C f {1+t- sr'i{\\h3{s)UinL- + ||V3/i3(s)||)ds, 
Jo 



\E\\<C{l + t)-^\\[uo,Eo,Bo\\L^nm 



+C [ il + t-s)-i\\h2is)\\L^nL^ds 
Jo 
ft 

72 J 



+ C [ {l+t~s) 
Jo 



■^i\\h3{s)Uir,L- + \\V'h3{s)\\)ds, 



and 



\\u\\<Cil + t)-i\\Uo\\L^nL- 

+C / {1 + t - s)--^\\hi{s)\\LinL^ds 
Jo 

+C f {1 + t - s)-^\[h2{s),h3(s)]UinL-ds. 
Jo 

It is straightforward to verify 

\\[hl,h2,h3\\LinL- + \\h3\\H^ < CSNiU). 

Then, ||i?|| is estimated by 

\\B\\<C{l+t)-i\\[uo,Eo,Bo\\L^nH^ 

+C [ {l+t-s)-i{l + syidsX{t) 
Jo 

< C{l+t)-l{\\[uo,Eo,Bo\\L^nH^ + X{t)), 
and in the same way, it holds that 

\\[u,E,B]\\<C{l+t)-H\\Uo\\L^nm+X{t)). 
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Therefore, one has 



f {l + .s)-i+'\\[u,E,B]fds 
Jo 



10 

< C Al + s)--^+^{l + s)-ids{\\U,\\l.^^.. + X{tf) 

JO 

<C{l + tn\\U,\\l,^H:.+X{tf). 

Substituting it into (|6.9p gives 

(6.10) (1 + t)i+'£N{U{t)) + / (1 + s)-^+'VN{U{s))ds 

Jo 

which imphes 

X{t)<C{\\Uo\\lrr,H- + r+X{tf). 

Since |lC/o||Lin_ff"+i is sufficiently small, X(t) is bounded uniformly in time and also 
(|6.7p holds true. This completes the proof Lemma l6Tl D 

6.3. Optimal large-time behavior. According to Corollarv lS.ll in the linear case, 
the time rate obtained in Lemma 16.11 is not optimal with respect to n, u and E 
in the solution U = [n,u, E, B]. In fact, the further estimates imply the following 
improved result. 

Theorem 6.2. Let N > A. Assume that initial data Uq '■— [uo^uq, Eq, Bq] satisfies 
()6.3p and ||[/o|lLin-ff"+^ ** sufficiently small. Let U = [n^u,E^B] be the solution to 
the Cauchy problem (|6.ip . (|6.2p . (|6.3p obtained in Theorem \6.1i Then, 

\\nm < {^ + tr'\\Uo\\L^nH-^^-, 
\Ht)\\ <{l + ty^\Uo\\L^nH^+-, 
\\Em <{l + t)-iH'S + t)\\Uo\\L^nH-^-, 
\\B{t)\\ <{l + t)-^\Uo\\L^nH-+^, 
for any t > 0. 

We shall prove the above theorem as follows. The estimate of B in fact has 
been obtained in Lemma |6. 11 the estimate of u will be given in Lemma 16.31 and the 
estimates of n and E will be given in Lemma 16.41 

Thus, we first consider the time decay estimate of u. Before that, we need a 
high-order energy inequality. Define 

£]^{U{t))^\\W[n,u,E,B]\\l^^^+Ki Y. (9"(7Vn),c>"M) 

l<|a|<Af-l 

+K2 Yl {d"VxE,d°'W xu) 

l<\a\<N-2 

+K1K2 Y (5"(-VxS),5"£;), 

2<|a|<W-2 

and 

T^wiUit)) = ||Vn||^„_i + ||V2u||^„_i + ||v2s||^«-3 + ||V^S||^„„4. 
Here, constants < ki, ^2 <C 1 are properly chosen as in (|2.12p . Then, we have 
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Lemma 6.2. //||C/o||h« is sufficiently small, then 

(6.11) j/NiUit)) + cV),{U{t)) < aWur + £UU{t)))£UU{t)) 

for any t > 0. 

Proof. Similar to obtain (|6.6p . from the system (|6.ip . one has 

p^{Uit))+cV%{Uit)) 

< Yl {{d"hud°'n) + {d°'h2, d°'u) + (a"/i3, d^'E)) 

l<\a\<N 

l<|Q|<Af-l l<|Q|<Af-l 

l<|Q|<Ar-2 l<|a|<Af-2 

-KiK2 X (a"VxB,9"/i3), 

2<|q|<JV-2 

which together with hi ~ ^g V • /13 imply 

(6.12) l£^^iU{t))+cV%{Uit)) 

< Y ((5"/ii,5"n) + (9"/i2,5"^i) + (9"/i3,5"ii;)) + !lV[/i2>3]ill,«-2. 

l<|Q|<Af 

It is easy to see 

mh2,h,]\\l.^,<£^{U{t)) [£^{U{t))+VUum . 

In what follows we estimate the right-hand first term of (|6.12p . Let 1 < |q;| < A^. 
First, it holds that 

{d^h3,d^E) = (9"(nw),5"i?) < ||Vn||H"-i|lVuiU.v-i|lV£;||H«-i. 
Next, for the term containing /12, one has 

(9"/i2,9"u) 

~ {d"{u ■ \/u),d°'u) + {d°'{nVn),d°'u) 

+ (9"(u X B),d°'u) + {d°'{nAu),d"u) 
< ||Vuj|^N-i j|V^uj|/f-N-i||Vuj|^iv-i + j|V?i||^jv-i||V^uj|/f-iv-i 

+ (||m|| • ||VS|| • ||V^w||iji + ||Vu||H"-i||V-B||^iv-i||V^w||/jjv-2) 

+ ||Vn||^iv-i||V2w||^«-i. 

Similarly, 

(a"/ii,5"n) -(a"V-(nM),a"n) < llVnll^^-iUV^Mll^iv-i. 
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By collecting the above estimates, it follows from (|6.12p that 

(6.13) ^^£UU{t))+cV%{U{t)) 

< (||^/,|| +fh (C/(t))V2)£h (t7(t))V2ph^(f/(i))l/2 

+ £]^(t/(t)) V2p5v(t/(t)) + £]^(t/(t)) [£]^([/(t)) + vuum . 

Since ||C/o||_f/" is sufficiently small and so is £iy(U{t)) uniformly for all i > by 
(|6.4p . (|6.1ip holds by applying the Cauchy inequality to (|6.13p . This completes the 
proof of Lemma 16.21 D 

Furthermore, wc define 

Y{t)= sup |(l + s)t[||u(s)|p+£]^([/(s))]|, t>0. 

0<s<t '- J 

Similar to obtain the uniform- in-time bound of X(t) in Lemma 16. 1[ we have the 
following result to show the boundcdness of Y{t) for all i > and thus the time 
decay rate of \\u\\. 

Lemma 6.3. // j|C/o|lLin_ff"+2 is sufficiently small, then 

(6.14) supr(f)<C||;7o||iinff«+- 

t>0 

Proof. The starting point is based on the high-order energy inequality ()6.1ip . Fix 
e > small enough. The time-weighted estimate on (|6.1ip gives 

(6.15) (1 + t)i+'£'},{U{t)) + / (1 + s)i+'V'%{U{s))ds 

JO 

<£n{Uo)+ f {l + s)i+^[\\u{s)f + £Uu{s))]£%{Uis))ds 
Jo 

+ / il + s)i+'£'^{U{s))ds. 
Jo 

It is easy to see that the right-hand second term is bounded by Y{t)'^ since 

/ (1 + s)3+^(l + s)-i-Us < C 
Jo 

for e > small enough. For the last term of (|6.15p . one has 

(1 + s)3+^f h ([/(s))d,s < / (1 + s)-^+'[VN+iiU{s)) + \\VBis)f]ds 

Jo 

< (1 + t)^||C/o|liinff«+^ + / (1 + s)i+'\\^B{s)\\^ds. 

Jo 

where ()6.10p was used. Here, from Corollarv l5.ll one can estimate ||Vi?j| as 
l|VS|l<(l + t)-^||(7o|UinH4 

+ / il + t-s)-i{\\h2is)\\LinL2 + \\Vh2is)\\)ds 
Jo 

+ I (l + t-.s)-^(||/i3(.s)||LinL^ + l|V%(.s)||)ds. 
Jo 
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Since 

\\h2{t)\\L^nL^ + \\^h2{t)\\ + ||/l3(i)||LinL= + l|V%(t)|| 

<\\Uo\\L^nH-^^il + t)-Hl+t)--sY{ty^^ 

<\\Uo\\L^nH-^^{l+tr'Yit)'/^, 
it follows 

l|Vi3||<(l + t)-ij|C/o|Uinff«+i(l + ni)'^')- 
Putting the above estimates into (|6.15p yields 

{l + t)i+'£%{U{t))+ f (l + s)t+^P5^([/(s))ds 

+ f (l + s)i+Hl + s)-i\\Uo\\lir,H-+^{l + Yis))ds 
Jo 

< (1 + tY\\UoWl^H-M^ + Yit)) + Y{t)^ 

which implies 

(6.16) sup (l + s)t4([/(,s))<||C/o||iinH"+^(l + ^W) +>'(*)'■ 

0<s<t 

The rest is to estimate the optimal time-decay of ||u|| in terms of Y{t). In fact, for 
||w||, from Corollary 15.11 one has 

\Ht)\\<{l+t)-i\\Uo\\Lir^L2+ [ il+t-s)-i\\hi{s)\\Lir,L-ds 

Jo 

+ / {l+t-syT\\[h2{s),h3{s)\\L^nL^ds 
Jo 

<{l+t)-^\Uo\\L^nH-+^Y{tY/^, 
that is, 

(6.17) sup (1 + s)i\\uis)f < \\Uo\\l^^H.+,Y(t). 

0<s<t 

Combining (|6.16p and (|6.17p gives 

Y{t) < C\\Uo\\l.^H-M^ + Y{t)) + Y{tf. 

Therefore (|6.14p holds true since ||C/o||LinH"+2 is sufficiently small. This completes 
the proof of Lemma 16.31 D 

To complete the proof of Theorem 16. 2| the rest is to consider the time decay 
estimates of n and E. In fact, we have 

Lemma 6.4. // ||C/oi|Lin_ff"+2 is sufficiently small, then 

(6.18) ||n(t)|| <C\\Uo\\L^nH-^'4l+t)-\ 

(6.19) \\E{t)\\ < C||C/oiUinH~+^(l + t)"^ ln(3 + t), 
for any t > 0. 
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Proof. For |jn||, from Corollary 15.11 one has 

Mm<{l + t)-i\\Uo\\L^nL^-+ [ il + t~s)-i\\[h2{s),h3{s)]U^nL^.ds. 

Jo 

Here, it holds that 

Then, (|6.18p follows. For ||£'(i)||, similarly, it suffices to notice from Corollary 15. II 
that 

Jo 

t 

{l + t- s)-t(|l/i3(s)|UinL^ + ||V2/i3(s)||)ds 



t 

-lj„llrr l|2 



<(l + t)-^||C/o||LinH^'+ / (l + t-s)-3(i + s)-M.s||C/o||iinH"+= 

Jq 

<l|C/o||LinH"+^(l + 0"^M3 + i)- 

Then (|6.19p holds. This completes the proof of Lemma W^ D 

Proof of Theorem ll.lt The global existence is proyed in Theorem 16. 1[ and the 

time decay rates arc obtained in Theorem 16.21 D 



7. Appendix 

In this section we proyide a deriyation of the Nayier-Stokes-Maxwell system p.ip 
from the one-species Vlasoy-Maxwell-Boltzmann system. The main tool that we 
shall use is the Liu- Yang- Yu's macro-micro decomposition, which was initiated in 
Liu-Yu |39] and deyelopcd in Liu- Yang- Yu [SQ , and also has extensiye applications 
in the study of the stability of the Boltzmann equation, cf. [6 5) and references 
therein. 

The Vlasoy-Maxwell-Boltzmann system is a kinetic model in plasma physics to 
describe the time eyolution of dilute charged particles (e.g., electrons and ions) with 
the self-consistent Lorentz forces [45]. Under the same physical assumptions as for 
the Nayier-Stokes-Maxwell system (|l.ip . the Vlasoy-Maxwell-Boltzmann system 
takes the form of 

{dtf + i-^J-{E + ixB)- Vj/ = Q(/, /), 



(7.1) 



Here, the unknowns are / = f{t,x,C) : (0,oo) x R^ x M^ -^ [0,oo), E = E{t,x) : 
(0, oo) X R3 ^ M3 g^^j 5 ^ ^(^f.^ 3,) . (Q^ o^-j X m3 ^ ^3^ ^^^Y^ j^^^ ^^ ^) standing for 

the number distribution function of one-species of particles (e.g., electrons) which 
haye position x and yelocity ^ = (^1)^2, ^3) at time i, and E{t,x) and B{t,x) 
denoting the electromagnetic field as in the fluid case. Notice that throughout this 



dtE-y,xB= / UdL 
dtB + y,y.E^ 0, 
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appendix, we still use ^, which is the same notation as in Section [3.21 to denote the 
micro velocity variable for simplicity, and this does not make any confusion. 

Let f(t, X, £,) be the solution to the system (|7.H) . As in [35], in terms of f{t, x, £,), 
we introduce five fluid quantities, that is, the mass density n(t,x), momentum 
density n(t, x)u(t, x) and energy density e{t, x) + \u(t, x)\'^ /2 by 



[pi 



n{t,x) 
n(t, x)ui{t, x) 
')]{t,x) 



h-\'' 



i 

MO.fit,^A)d^, * = 1,2,3, 



and further define the local Maxwellian as 

_ n{t,x) 



M = Mr. 



,em 



{2Tre{t,x))3/' 



■ exp 



\^-uit,x)\-^ 

2e{t,x) 



Here 9{t,x) is the temperature related to the internal energy e{t,x) by e = |6', 
u{t,x) = (ui(t,x), U2(t, a;), U3(t,a;)) is the fluid velocity. As usual, tpai^), ex. — 
0,1,2,3,4, are the five collision invariants: 

^,(e) = e., « = 1,2,3, 

With the above definitions, the Maxwell equations in system (|7.ip can be rewritten 
as 



(7.2) 



dtE — Vx X i? = nu, 

dtB + V^xE^ 0, 

\/ ■ E ^ n^ - 71, Vx- B ^Q. 



This is exactly the same as the part of the Maxwell equations in system (|l.ip . The 
rest is to derive the same as the part of the Navier-Stokes equations in system (|l.ip . 
As in [3S], let us write f{t,x,£,) as the sum of the local Maxwellian M and 
microscopic component G = G{t,x,^), that is. 



(7.3) 



/(t,x,0=M(t,x,0 + G(t,x,0. 



By doing this, the system (j7.1|) can be reformulated into a coupled system containing 
the conservation laws for the macroscopic fluid components and the equation for 
the microscopic component. We shall sec that the Navier-Stokes-Maxwell system 
(|l.ip turns out to be the macroscopic fluid part. 

For later use, we need to define the projection with respect to the Maxwellian 
M. To do that, define the inner product 



)m over L^ as 



{h,9) 



M 



h{Og{OM-'dt 
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As in [m, the macroscopic linear subspace M spanned by {-00(0^; V'i(OM, i 
1, 2, 3, V'4(C)M} has an orthogonal basis consisting of 



Xo 






^^^W^' *-l'2'3, 



X4 = :^^^~3 M, 



, {Xa.Xp) 



M 



5^B: fora, /?-0,l,2,3,4. 



The macroscopic projection Pq and microscopic projection Pi can then be defined 

by 

4 

Voh = ^{h,Xa)^Xa, Pih = h-Poh. 



Q=0 



Now, we go back to the Vlasov-Maxwell-Boltzmann equation (|7.ip . The property 
of collision invariants implies 

for a = 0, 1, • • • ,4. Using the decomposition (|7.3p fm-ther gives 
' Stp + V^E • (nu) = 0, 



(7.4) 



dt{nui) + X) dj{nuiUj) + a^P + n(_B + m x _B)j 

+ / V'.(O^V,Gde = 0, * = 1,2,3, 
P (^ + s)] + E d] {% [P (^ + e) +p] } + "w • (S + u X B) 

+ / MO^ ■ v.Gdc = 0. 

Here and hereafter we use di = dx^ for simplicity. In (|7.4p . the state equation for 
the monatomic gas is given by P = |ne and hence P = nO. By noticing that 9tG, 
V^G, LmG, and Q{G,G) are microscopic, the microscopic equation is obtained 
by applying the microscopic projection Pi to the first equation of system (|7.ip : 



(7.5) dtG + Pi(C • V,M + C • V,G) -{E + ixB)- V^G = LmG + 0(G, G), 

where Lm is the linearized collision operator defined by 

LM.g = Q(M + 5, M + .g) - 0(.g, .g) = Q(M, g) + Q{g, M). 

It is known that the null space of Lm is exactly M. Since Lm is a bounded and 
one-to-one operator on M-'- , from (|7.5p , one has 



(7.6) 
with 



G = Lj^i (Pi (C • V,M)) + e 



(7.7) 



e = Lj^' {dtG + Pi (e • V,G) -{E + ixB)- V^G - Q{G, G)) 
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Substituting (|7.6|) into ()7.4|) yields the following fluid-type system 
' dtn + Vx ■ {nu) = 0, 



(7.; 



dt{nui) + J2 ^j [nUiUj) + S^P + n(i? + u x i?)^ 

= E d, He)n,] - / 7/;,(e)e • v.ede, ^ = 1,2,3, 

9t [n(i|u|2 + e)] + V^ • {u{n{\\u\'^ + e) + p)) + tiu • (S + u x B) 

1,7 = 1 7 = 1 JM3 



with 



Tij = djUi + diUj - -5ij\/ X ■ u- 



Here, the viscosity coefficient v{9) and heat conductivity coefficient k{9) both de- 
pending only on 9 can be explicitly given by introducing the Burnett functions, 

cf. Eg. 

Observe that the microscopic term Q in (|7.7p is generated by the linear super- 
position of the first-order time-space derivative of G and the nonlinear terms. If 
all the terms containing 8 in (|7.8p are dropped, then by recalling (|7.2p . the fluid 



quantities n, u, e = |0 and the electromagnetic field E, B satisfy the following 
coupled Navier-Stokes-Maxwell system in the non-isentropic case 

' dtn + \/x ■ {nu) ~ 0, 

3 

dtinui) + Y^ dj (nUiUj) + diP + n{E + u x B)i 

= i:dj[y{e)n.j\, i- 1,2,3, 
dt [n(i|u|2 + e)] + V^ • {u{n{\\u\'^ + e) + p)) + nu ■ {E + u x B) 

= E d, M9)u.nj} + E d, {n{9)dj9) , 
4,i=l i=l 

dtE — Vx X B = nu, 

dtB + \'xXE = 0, 

V • £■ == rib - n, \7x-B = 0. 

System (|l.ip that we have studied here is the simplified form in the iscntropic case. 
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